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Preface 


This report is the User’s Manual for the Metal Matrix Laminate Tailoring (MMLT) 
program. The code is capable of tailoring the fabrication process, constituent characteristics, 
and laminate parameters (individually or concurrently) for a wide variety of metal matrix 
composite (MMC) materials, to improve the performance and identify trends or behavior of 
MMCs under different thermo-mechanical loading conditions. This document is meant to 
serve as a guide in the use of the MMLT code. Detailed explanations of the composite 
mechanics and tailoring analysis are beyond the scope of this document, and may be found in 
the references. MMLT has been developed by the Structural Mechanics Branch at NASA 
Lewis Research Center (LeRC). 

Please note that the MMLT code is being made available strictly as a research tool. Neither 
NASA LeRC, nor any contractors, nor grantees that have contributed to the code 
development, assume liability for application of the code beyond research needs. 

Any questions or related items concerning this computer code can be directed to Dr. P. L. 

N. Murthy at the NASA LeRC, Cleveland, OH 44135 or call (216) 433-3332. 
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INTRODUCTION 


1.1 Scope 


Future advancements in the propulsion and aerospace industry require a dramatic extension in 
materials and processing technology. Current research and development suggest that metal 
matrix composites (MMCs) will be essential in meeting performance goals in civil transport 
systems that require increased fuel economy, improved reliability, extended life, and reduced 
operating costs. However, there is a demand to generate the technology needed that will 
enable the development of MMCs to meet these lofty goals. MMLT (Metal Matrix 
Laminate Tailoring) is a Fortran computer code, developed at NASA Lewis Research Center 
to aid the researcher in producing some of the technology required to advance the 
fabrication, selection, and tailoring of MMCs as functional materials in the design of high 
temperature engine components. 

Through computational simulation of the current fabrication process and thermo- 
mechanical (TM) loading, the MMLT code is capable of tailoring the fabrication process, 
material characteristics, and laminate parameters, individually or simultaneously, for MMCs. 
By tailoring these parameters the user has the ability to: control the microstresses during 
fabrication [Reference 1] and subsequent TM loading [Reference 2]; control the development 
of mechanical properties and micromechanical parameters at the constituent, ply, and 
laminate level; and improve the TM load carrying capacity of the fabricated composite 
[Reference 1-2]. 

The basis of the MMLT code encompasses two computer codes, METCAN (METal 
matrix Composite ANalyzer) [Reference 3] and the ADS optimizer [Reference 4-5]. To 
capture the nonlinear behavior of a MMC during fabrication and the subsequent TM loading, 
nonlinear composite mechanics are used as encoded in METCAN. The composite mechanics 
include effects of the elevated temperature, the nonlinear response of the constituent 
materials, and the residual stress build-up during fabrication and TM loading. In order to 
handle the multifarious optimization problem, the tailoring of the MMC is formulated as a 
nonlinear programming problem and is solved with the feasible directions method. By taking 
advantage of the unique capabilities of these two methods, the foundation of the MMLT code 
has been formed. 
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1.1 Scope 


MMLT provides the unique feature to either simulate or tailor the complex behavior 
of MMCs from its fabrication process through its in-service life by considering several 
design parameters simultaneously. The potential user will be able to determine trends for 
fabrication and composite parameters before expensive experimental investigations are begun. 
Also, further insight of the material system under different loading conditions can be 
quantified. Results from the MMLT package can be used in the design and analysis of 
composite structures by tailoring both the fabrication process, constituent parameters, and 
laminate characteristics for specific requirements. However, the potential user of MMLT is 
reminded that the program is in an evolutionary state and the methodology which MMLT 
comprises is of an on-going nature. 

The remaining sections of Part I include a brief background description of MMLT 
(which includes most of the assumptions and limitations pertaining to this code), an example 
input file with instructions, presentation of the geometric conventions used, and input/output 
files employed and produced by MMLT. Part n contains a detailed description of the 
command names needed to create a data file that can be executed by the MMLT code. Also, 
Part II has two example data files to aid the user in creating their own input files. Part HI 
exp lain s the procedures required to execute the MMLT code after the data file has been 
created on either the VAX, Cray X-MP, or Convex computer systems which are available at 
NASA Lewis Research Center. 
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1 .2 MMLT Overview 


In order to understand the motivation behind the MMLT code, a brief background is 
presented so that the user can apply the methodology to enhance the performance of the 
material. The objective of this section, therefore, is to identify existing problems within 
MMCs, present the physical model used in the simulation, and describe the methodology 
used to simulate and tailor the complex behavior of MMCs. 

Applications of MMCs are becoming a more viable solution for high material 
property demands in terms of high temperature operation, high modulus, strength, hardness, 
conductivity, dimensional stability, and low density, but MMCs still have not reached their 
full potential because of a lack of technology. Furthermore, extensive experimental 
investigations are prohibitive in time and cost, as a result, it is advantageous to have 
computational schemes which can simulate the behavior of MMCs as a viable alternative to 
rigorous experimental investigations. This will also minimize the number of experiments to 
be performed by eliminating unwanted candidate MMCs. The MMLT code was developed 
to tailor the fabrication process, interphase layer characteristics, and laminate parameters, 
either individually or concurrently, in order to improve the thermo-mechanical performance 
and extended the service life of MMCs. The code can be used to develop trends and lend 
insight into the impact of different tailorable parameters for MMCs so that the overall 
performance can be improved. 

Significant problems exist within MMCs, for example, the presence of residual 
stresses in the composite and the interfacial bonding conditions between the fibers and 
matrix. The residual stresses are developed during cool-down in the fabrication phase 
because of the mismatch in the coefficients of thermal expansion (CTE) between the fiber 
and matrix. In composite laminates, additional residual stresses are developed because of 
differences in the CTE between plies. Depending on the properties of constituents, the fiber 
volume ratio, and the temperature differential involved, the residual stresses may reach 
critical magnitudes which degrade the mechanical properties and thermo-mechanical (TM) 
load carrying capacity of the composite. This may even cause local damage, such as matrix 
cracking, fiber debonding, and degradation of the constituent properties. Therefore, the 
development of residual stresses is crucial to the performance of many MMCs. 

Although these critical residual stresses exist within MMCs, it is very possible to 
control them by: (1) altering critical fabrication parameters, such as processing temperature 
and pressure histories; (2) incorporating a compatible interphase layer between the fibers and 
the matrix; (3) modifying laminate parameters, such as fiber volume ratio and ply 
orientation. Candidate design variables include temperature and pressure histories of the 
fabrication process; the characteristics of a candidate interphase layer, such as, the modulus, 
CTE, strength, and layer thickness; and laminate parameters, such as the fiber volume ratio, 
ply thickness, and ply angles. The user is able to control the build-up of microstresses 
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1.2 MMLT Overview 


during fabrication and subsequent TM loading through stress constraints, while ensuring the 
integrity of the tailored MMC by avoiding failure at the constituent, ply, and laminate levels. 
Also, improvements in the mechanical properties of the tailored composite and an increase in 
the load carrying capacity of a MMC system can be achieved. 

A typical thermo-mechanical life cycle of a MMC laminate is schematically shown in 
Figure 1.2-1. The laminate is fabricated consolidating the raw materials at an elevated 
temperature followed by gradual cooling to room conditions. Consolidation pressure is 
applied to diffuse and bond the matrix with the fibers. The combination of temperature and 
pressure during fabrication results in a hot isostatic pressing (HIP) process. The residual 
microstresses are developed during the cool-down consolidation phase of the fabrication 
process. 

The composite behavior during fabrication and TM loading is computationally 
simulated with incremental nonlinear composite mechanics [Reference 3]. The composite 
mechanics include inelastic, thermal and residual stress buildup effects on the constituent, 
lamina, and laminate levels. In this manner, the procedure provides more accurate 
predictions and is able to capture the coupling between the fabrication process, interphase 
layer characteristics, and laminate parameters. The mechanics incorporate a multi-cell model 
packed in a square array. An idealization of a single unit cell is shown in Figure 1.2-2 with 
the three material phases available; fiber (f), matrix (m), and an optional interphase (d). 

For the optimization, the ADS optimizer [Reference 4-5] is used. A constrained 
nonlinear programming (feasible directions) method is utilized to minimize an objective 
function, F(Z). The problem solved, in the most general form, is: 

Minimize F(Z) 

Subject to constraints of the following form: 

Z L < z < z u 
Q(Z) < 0 

The design variables are represented by the vector Z, with U and L indicating the upper and 
lower bounds. Also, Q(Z) are inequality performance constraints. The feasible domain is 
defined by the last two equations with a one dimensional search being performed. Scaling 
and normalization of the design variables, objective function, and constraints are all done 
internally in order to improve the efficiency, reliability, and the ease-ability to the user. For 
further information on the ADS optimizer the user is referred to the references [4-5]. 

Displayed in Figure 1 .2-3 is the modular format of the MMLT code. A user friendly 
input format is provided to the user. After input is read in, a readable METCAN file is 
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1.2 MMLT Overview 


created from the MMLT. In turn, the analysis of the composite is performed using 
composite micromechanics and laminate analysis. The results from the METC AN module 
are used to calculate the material performance measures needed by the ADS optimizer 
module. Upon completion of the optimization, output files are provided to the user. 


Figure 1.2-1: Typical Fabrication and thermo-mechanical loading of metal matrix laminates 
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1.2 MMLT Overview 


Figure 1.2-2: Idealization of metal matrix composite unit cell 



Regions of Constituent 
Material Nonuniformity 

1 


A: Matrix 



B: Matrix and Interphase 
C: Matrix, Fiber, and Interphase 
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1.2 MMLT Overview 


Figure 1.2-3: Modular Structure of MMLT 



Metal Matrix Laminate Tailoring 
NASA Lewis Research Center 


Part / - Introduction 

1.2-5 








1 .3 MMLT Input Example 


A sample case with a step by step description of the procedure (Table 1.3-1) required to 
prepare MMLT input data file (Table 1.3-2) is provided in this section. The objective of this 
example case is to tailor the fabrication process for a unidirectional graphite (P100)/copper 
(Cu) MMC in order to reduce the longitudinal and transverse matrix microstresses at the end 
of fabrication. These high stresses may cause degradation in the composite properties, fiber 
debonding, and reduce the thermo-mechanical fatigue life. Therefore, it may be beneficial to 
control these stresses during fabrication. 

Current temperature and pressure histories are shown in Figures 1.3-la and 1.3-lb. 

A total of 5 load history points (i.e., the control points defining the temperature and pressure 
history) are needed to simulate the fabrication process and the five points are marked with 
their appropriate coordinates. As a result, four entries are needed for the subincrements (the 
number of steps between the load history points), eg. from point one to point two is the first 
entry of the subincrement, from point two to point three a second entry of the subincrement 
is needed and so forth. Each designated point on the plots is a design variable except for the 
consolidation and room temperatures (the first and final temperatures) and the final pressure 
which is set to zero. Since the build-up of residual stresses and the integrity of the 
composite material are primarily affected by the temperature and consolidation pressure 
histories only the cool-down phase of the fabrication process is simulated [Reference 1]. 

The fiber volume ratio (FVR) is 40%, no voids are present, the ply thickness is 0.04 
in., and a perfect bond (no interphase) is assumed between the fiber and matrix. Stress 
constraints (the current stress is greater than the current compressive strength and less than 
the tensile strength at any time during the fabrication) on the fiber and matrix are activated 
throughout fabrication. This ensures the integrity of the tailored MMC by avoiding failures 
at the micro-level. 

The individual steps, along with the corresponding line numbers in the example input 
data file (refer to Table 1.3-2), that are required to model this problem are presented in 
Table 1.3-1. The temperature and pressure load histories are taken from Figures 1.3-la and 
1.3-2b. Line numbers are provided in the sample input deck for reference purposes only and 
are not part of the input file. 
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1.3 MMLT Input Example 



Step No. 

Line No. 

1 

1 

2 

2 

3 

3 



Des cription of Data 

Use a comment, record to describe the input file, optional. 

Use the 'PROCESS' record to indicate the load history points (5 

control points), which define the temperature and pressure histories in 
Figures 1.3-1 a and 1 .3-1 b. Also the number of steps (subincrements) 
between the linear segments of the load history. The number of 
entries of subincrements is one less than the number of control points 
defining the load history. For example, from control point one to 
control point two the first entry for the subincrement is six, i.e. there 
are six steps between the first linear segment designated by control 
points one and two. 

Use the 'TIME' record to specify the time range at each control point 
over the entire process. A numerical label is placed in the 4th field 
preceded by the mnemonic record and two blank fields. A total of 5 
points are needed to match the load history points in the PROCESS 
record. 

Use the 'TEMPR' record to specify the temperature history during the 
cool-down phase of fabrication. Refer to Figure 1 .3-1 a for the current 
values at each control point. 

Use the 'PRESR' record to specify the pressure history during the 
cool-down phase of fabrication. Refer to Figure 1 .3-1 b for the current 
values at th e corresponding control points. 

Use the 'NY' record to specify the transverse applied force at each 
load segment. The transverse force is equivalent to the corresponding 
pressure value (refer to 'PRESR' record) multiplied by the laminate 
thickness (refer to the 'MATERIAL' record). This is used to simulate 
the ho t isostatic pressing of the fabrication process. 

The 'LAY-UP' record defines the order of the ply lay-up by a user 
supplied name. In this example the composite is unidirectional, thus, 
onl y one ply (PLY1) is defined. _____ 

Use the 'MATERIAL' record to define the material of the ply matching 
the name (PLY1 ) given in the 'LAY-UP' record. Additional data in this 
record are the fiber and matrix code names used to retrieve the 
constituent material properties from the databank. For the graphite 
fiber use PI 00 and COPR for the copper matrix, the FVR (0.40), the 
interphase thickness (0.0) as a percent of the fiber diameter, the ply 
thickness (0.04), and the fiber orientation angle (0.0). 
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1.3 MMLT Input Example 


| Step No. 

Line No. 

Description of Data 

9 

9-15 

Use the 'VARIABLE' record to define ail design variables individually 
by assigning a user supplied name (e.g., TEMP2), the corresponding 
variable number (e.g., 2 for variable TEMP3) in increasing order, the 
position of variable in the global design vector (refer to Table 2.2-1 in 
Section 2.2.3 pg. 2.2-21), the lower bound, the upper bound, and the 
initial value of the design variable. For this example there are seven 
design variables, three temperature and four pressure values. 

10 

16-19 

Use the 'DEPEND' record to define the four dependent design 
parameters to the pressure design variables defined in the previous 
record. This defines the consolidation pressure in the y direction 
required by the fabrication process to consolidate the composite. A 
user supplied name is given for each parameter and the position in the 
storage vector is given from Table 2.2-1 in Section 2.2.3. The 
storage position of the design variable that the dependent design 
parameter is linked to, is also given. Finally, a ratio to the design 
variable is defined. 

11 

20-37 

Use the 'CONSTRNT' record to define the constraints by assigning the 
user supplied name (e.g. SF11C), multiplication factor (1.0), equation 
form of the constraint (e.g. 1, refer to Table 2.2-2 in Section 2.2.3 
pg. 2.2-24), type of constraint (e.g. 0, nonlinear inequality constraint), 
and location of constraint in the storage array, two entries are used to 
specify the beginning and end point within the storage array, refer to 
Table 2.2-3 in Section 2.2.3 pg. 2.2-25. A second record is needed, 
the '&CNST' mnemonic, which is an extension of the first record. 

This second record contains the bound needed to evaluate the 
constraint equation chosen by the user from Table 2.2-2. 

12 

38 

Use the 'OPTMIZE' record to activate the optimization solution by 
defining the optimizer strategy (0), the optimization procedure (5), the 
search procedure (8), the output request (3440), the scaling indicator 
(1), the number of design moves (40), and the maximum number of 
optimization evaluations (200). Except for the last two values, all 
other values are default values and should be used. Consult 
Reference 5 for further details. 

13 

39 

Use the 'OBJTYPE' record to specify the objective function. For this 
example, the longitudinal and transverse matrix microstresses are 
minimized, therefore, 3 is chosen as the objective function (refer to 
pg. 2.2-29 for a list of available objective functions). 
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1.3 MML T Input Example 


Step No. 


Line No. 


Description of Data 


14 


40-41 


Use the 'SENSE' record to indicate the stresses that are to be 
minimized in the objective function. Since both matrix microstresses 
are considered in absolute values and weighted evenly, a positive 1 .0 
is inserted in the proper place to activate o m11A and o m22A with 0.0s in 
the remaining fields to nullify the other matrix microstresses. A 
continuation ' + ' record is used, though not all matrix microstresses (a 
total of nine) are used in the objective function. An optional 
continuation card is used for completeness. 


15 


42 


Use the 


'*END_INP' record to indicate the end of input. 








1.3 MMLT Input Example 


Table 1 .3-2: MMLT Sample Input 

Line No. 1 2 3- 4 5 6 7 8 9 

1 $ Tailor the temperature and pressure histories of a [0] PlOO/Copper MMC. 


2 

PROCESS 

5 

6 

35 

10 5 




3 

TIME 



1 

0.0 1800.0 

3600.0 

5400.0 

5700.0 

4 

TEMPR 



2 

1500.0 1500.0 

1075.0 

650.0 

70.0 

5 

PRESR 



3 

0.0 -1000.0 

-1000.0 

-1000.0 

0.0 

6 

NY 



4 

0.0 -40.0 

-40.0 

-40.0 

0.0 

7 

LAY-UP 

PLY1 







8 

MATERIAL 

PLY1 

P100COPR 

0.40 

^3* 

O 

• 

o 

o 

• 

o 

0.0 



9 

VARIABLE 

TEMP 2 

1 

12 

70.0 

1500.0 

1500.0 


10 

VARIABLE 

TEMP 3 

2 

13 

70.0 

1500.0 

1075.0 


11 

VARIABLE 

TEMP4 

3 

14 

70.0 

1500.0 

650.0 


12 

VARIABLE 

PRES1 

4 

HI 

-25000.0 

0.0 

0.0 


13 

VARIABLE 

PRES2 

5 

112 

-25000.0 

0.0 

-1000.0 


14 

VARIABLE 

PRES3 

6 

113 

-25000.0 

0.0 

-1000.0 


15 

VARIABLE 

PRES4 

7 

114 

-25000.0 

0.0 

-1000.0 


16 

DEPEND 

NY1 

201 

111 

0.04 




17 

DEPEND 

NY2 

202 

112 

0.04 




18 

DEPEND 

NY 3 

203 

113 

0.04 




19 

DEPEND 

NY 4 

204 

114 

0.04 




20 

CONSTRNT 

Sf 22C 

1.0 

1 

0 51 

55 



21 

&CNST 

0.85 







22 

CONSTRNT 

Sm22A 

1.0 

1 

0 151 

155 



23 

&CNST 

0.85 







24 

CONSTRNT 

Sm22C 

1 « 0 

1 

0 351 

355 



25 

&CNST 

0.85 







26 

CONSTRNT 

SmllA 

1.0 

1 

0 451 

455 



27 

&CNST 

0.85 
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1.3 MML T Input Example 


Line No. 

1- 

2 

3 

4- 

28 

CONSTRNT 

SmllC 

1.0 

1 

29 

&CNST 

0.85 



30 

CONSTRNT 

SfllC 

1.0 

1 

31 

&CNST 

0.85 



32 

CONSTRNT 

all 

1.0 

2 

33 

&CNST 

0.85 



34 

CONSTRNT 

a22 

1.0 

2 

35 

ScCNST 

0.85 



36 

CONSTRNT 

a33 

1.0 

2 

37 

ScCNST 

0.85 



38 

OPTIMIZE 

0 

5 

8 

39 

OBJTYPE 

3 



40 

SENSE 

1.0 

1.0 

0.0 

41 

+ 

0.0 



42 

*END INP 





5 6 7 8 9 

0 651 655 

0 751 755 

0 1931 

0 1932 

0 1933 

3440 1 40 200 

0.0 0.0 0.0 0.0 0 . 0 + 
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Temperature (F) 


1.3 MMLT Input Example 


Figure 1.3-1: Current fabrication process for a PlOO/Cu Unidirectional composite 



Time (sec) 

(a) 



(b) 
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1 .4 MMLT Flowchart 


The computational procedure for the solution of the tailoring methodology is schematically shown in 
Figure 1.4-1. The input file is read to determine the fabrication process, the thermal-mechanical 
loads applied after fabrication, constituent properties, and ply data. Also, optimization parameters 
are needed, such as design variables, constraints, and designation of the objective function. The 
MMC model is then analyzed by METCAN to determine the response of the composite at the 
constituent, ply, and laminate level under the loading history. Examples of the response are the 
micro/macro stress histories throughout fabrication and TM loading, and material properties at the 
constituent, lamina, and laminate level. 

The response of the composite material, during fabrication and/or after the TM load is 
applied, is used to calculate the objective function and constraints chosen by the user at the time of 
input. By considering the sensitivity of constraints and the size and location within the feasible 
domain, the optimizer will determine new design variables to reduce the objective function while 
remaining within the acceptable constraint limits and feasible region. Because of the large number 
of design variables and constraints used in MMLT and the nonlinearity of the problem, the method 
of feasible directions is used. The concept behind the method of feasible directions is to rapidly 
reduce the objective function while maintaining a feasible design. The sensitivity of the design (i.e. 
the improvement of the design variables) is linked to the objective function and constraints chosen 
by the user at input. If convergence between the current and new design variables and objective 
function are within an acceptable tolerance, the program will cease operation, if not, another 
iter ati on will be performed till an acceptable design is found. For further details on the composite 
response or the ADS optimizer used in MMLT, the user should refer to references [3-5]. 
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1.4 MMLT Flowchart 


Figure 1.4-1: Flow Chart of MMLT 
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1 .5 MMLT Geometric Conventions 


The geometric conventions used in MMLT are similar to conventions used in the METCAN code 
[Reference 3], but for completeness will be repeated herein. To model the composite material, the 
global and local coordinate systems and the conventions for ply ordering, measuring orientation 
angle, and loading are presented. The English Engineering System is used throughout this manual 
and the example problems, but the SI units can be used if all entries are consistent, i.e. in the 
primary MMLT input file and the constituent databank file. Refer to Appendix A for the 
nomenclature and unit values used in this manual. 


1.5.1 Global (Laminate) Coordinate System 

For an arbitrary laminate, shown in Figure 1.5-1, the global coordinate system is defined. An 
unidirectional laminate has its fibers running parallel to the x-axis, whereas the y-axis is orientated 
perpendicular to this fiber direction, and the z-axis is orientated through the thickness of the 
composite. 


Figure 1.5-1: MMLT Global (laminate) coordinate system 


Z 
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1.5 MMLT Geometric Conventions 


1.5.2 Local (Material) Coordinate System 

The local (material) coordinate system on a composite ply is shown in Figure 1.5-2. The direction 
parallel to fibers is called the longitudinal direction (axis 1). The direction perpendicular to the 
fibers is called the transverse direction (axis 2). The direction for thickness of the ply is the 3-axis. 


Figure 1.5-2: MMLT local (material) coordinate system 
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1.5 MMLT Geometric Conventions 


1.5.3 Ply Ordering Convention 

The ordering of the plies is depicted in Figure 1.5-3. The first ply is placed at the bottom or in 
plane z=0. Subsequent plies are stacked on top of the first ply in the positive z direction. 


Figure 1.5-3: MMLT ply ordering convention 
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1.5 MMLT Geometric Conventions 


1.5.4 Ply Angle Orientation 

The angle $ of a ply is the angle between the local and global coordinate systems. 6 is taken positive 
when the angle of the 1-2 axes measured from the x-y axes is in the counterclockwise direction as 
shown in Figure 1 . 5 - 4 . 


Figure 1.5-4: MMLT ply orientation angle convention 
2 



X, Y: global coordinate axes 

1 , 2: local coordinate axes 

0: angle between global and local axes 
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1.5 MMLT Geometric Conventions 


1 .5.5 Loading Conventions 

The various mechanical loadings (pressure P, resultant force N, and resultant moment M) and then- 
positive sense are presented in Figure 1.5-5. Units for the pressure on the upper and lower surface 
are force per area, where N„ N y , and have the units of force per unit depth, and M x , M y , and 
Mjy the units of moment per unit depth. All loading is applied at the global level. 


Figure 1.5-5: MMLT positive sense of pressure, resultant forces, and resultant moments 


Mxy 
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1 .6 MMLT Input/Output File Structure 


The file structure of the input/output for MMLT is shown in Figure 1.6-1. Two input files are 
needed to execute the MMLT code, a primary input file and a constituent databank file. MMLT 
produces four output files: (1) a standard METCAN output file; (2) an input file to METCAN; (3) 
ten postprocessing files; and (4) a history of optimization. If further details are required on the 
constituent databank file or any METCAN output file, refer to Reference 3. 

The primary input file contains all the fabrication processing history, thermal-mechanical 
loading, and identifies the constituent materials and the optimizer data. Data needed by the 
optimizer includes design variables, constraints, optimization strategy, and objective function. The 
second file is a databank of constituent material properties needed to begin the analysis. By 
specifying the constituents in the primary input file, their associated reference properties are found 
in the databank. A description of the constituent databank is located in Appendix B. 

Output files include a standard METCAN output file, a METCAN input file, and ten 
METCAN postprocessing files related to optimum designs. The output files generated by the 
METCAN module describe the simulation of the tailored MMC and comprises information at three 
levels, the laminate, ply, and constituent levels. Figure 1.6-2 shows the three levels of analysis and 
the various information associated with each level. A METCAN input file is provided to the user if 
future analysis is needed on the tailored MMC. Figure 1.6-3 provides an outline of the METCAN 
primary input file and a listing of information contained in the file and additional information that 
can be added. If MMLT was executed for one step, the MMLT input deck would be equivalent to 
the METCAN input file. 

The ten postprocessing output files are automatically generated by MMLT. These ten files 
essentially contain the same information generated in the primary METCAN output file, but in a 
more structured format for the user’s convenience. The structured format of the files can be very 
easily manipulated for postprocessing the results, e.g., generating a stress/strain curve for the 
tailored MMC. Figure 1.6-4 describes the ten postprocessing output files and their contents. 

A history of the tailoring steps and a detailed description of the optimization search is located 
in an output file containing all the steps taken by the optimizer. For example, all design variables 
are displayed at each step, the active and violated constraints, loading history, a select number of 
constituent microstresses at the end of fabrication and TM loading, and material properties. The 
output file from the optimizer is a summary of important information needed by the user to evaluate 
the results. A summary of the results produced by the optimizer module is shown in Figure 1.6-5. 
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1.6 MMLT Input/Output File Structure 


Figure 1.6-1: MMLT Input/Output File Structure 
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1.6 MMLT Input/Output Fite Structure 


Figure 1.6-2: Levels of Analysis Information Available in the Primary Output File. 


Primary output file 

Output reporting tailored by the user 


Composite 

• 2-D and 3-D material 
properties 

• 3-D stress-strain relation 

• Force-displacement 
relations 

• Reduced membrane and 
bending stiffnesses 

• Stress concentration 
factors around a circular 
hole in an infinite plate 

• Elastic and material 
properties for a finite 
element analysis 


P!y 

• Load/time step history 

• Failure index for each ply 

• Ply properties 

• Ply stress-strain relations 
in local coordinates 


Constituents 

(fiber, matrix, and interphase) 

• Echo reference 
constituent material 
properties found in the 
databank 

• Instantaneous values of 
constituent properties 

• Microstresses and micro- 
strains for each 
constituent in local 
coordinates 
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1.6 MMLT Input/Output File Structure 


Figure 1.6-3: Functions of the primary input file 
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1.6 MMLT Input/Output Fite Structure 


Ten Postprocessing Files 
Contains the same information 
presented in METCAN output 
file, but formatted for user 
designed postprocessing. 


| Constituent Level 

Ply Level 

Composite Level 

File 2 

Time and Instantaneous 
Fiber Properties 

File 1 

Time and Temperature 
for Each Ply 

File 10 

Time and Instantaneous 
Composite Properties 

File 3 

Time and Instantaneous 
Matrix Properties 

File 5 

Time and Instantaneous 
Ply Properties 


File 4 

Time and Instantaneous 
Interphase Properties 

File 9 

Time, Ply Stresses, Ply Strains, Composite Stresses, 
and Composite Strains 

File 6 

Time, Fiber Microstresses, 
Fiber Microstrains 

File 7 

Time, Matrix Microstresses, 
Matrix Microstrains 


File $ 

Time, Interphase 
Microstresses, 
Interphase Microstrains 


Figure 1.6-4: Contents of the Ten Postprocessing Files. 
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1.6 MMLT Input/Output File Structure 


Figure 1.6-5: List of information available in the history of tailoring output files 
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MMLT INPUT PREPARATION 


2.1 General Description 


Input to the MMLT program is provided in a straightforward user interface. The ASCII input file is 
created by using "record groups" identified by individual character mnenomics followed by the 
pertinent information. Character mnemonics are organized in a specific order. Therefore, each 
record in the input file should be defined in the order they are found in the following sections. The 
Hata deck is divided into four different sections. The organization of input mnemonics for MMLT is 
shown in Figure 2.1-1. 

The commands for MMLT are divided into four sections of input. Section 2.2. 1 contains the 
information needed to simulate the fabrication process of the metal matrix composite and thermal- 
mechanical loading, i.e., the time, temperature, pressure, and mechanical loading at each load step. 
In the next section (2.2.2), the composite material identification and laminate parameters are 
defined. Section 2.2.3 assigns the design variables and constraints needed for the optimizer module. 
The subsequent section (2.2.4) defines the optimization parameters and the objective function. 
Typical examples of MMLT data files are located in Section 2.3. 
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2. 1 General Description 


Figure 2.1-1: Organization of Input File for MMLT 


PROCESS 

TIME 

TEMPR 

PRESR 

NX 

NY 

NXY 

MX 

MY 

MXY 


PRIORITY 

LAYUP 

MATERIAL 


Fabrication 

and 

Thermal-Mechanical Loading 


Laminate Lay-up 
and 

Material Identification 


VARIABLE 

CONSTRNT 


Design Varialbes and Constraints 


OPTIMIZE 

OBJTYPE 

SENSE 


Optimization Parameters, 

Objective Function, 
and Weighing Factors 


*END_INP 


End of Input 
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2.2 MMLT Data Deck 


This section describes the identification mnemonics needed to create the data deck for MMLT. The 
command name is at the top of the page followed by a description. The proper format and an example 
is also included. All fields (ten per line) are eight spaces in length with the first field being a character 
field (A8). Ensuing fields can be character (A8), real (F8 or E8) or integer (18) and are specified in 
the data deck sections. Alphanumeric and integer data can be left or right justified, while real data can 
be entered anywhere in the appropriate field. Two general commands are given which can be used 
anywhere within the data deck, i.e., the comment and continuation records. 
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2.2 MMLT Data Deck 


$ 


Description: 

The $ symbol is used to make comments within the data deck. Also, the $ can be used to inactivate 
the current line. 


Format and Example: 

$ COMMENT 

1 2 3 4 5 6 7 8 9 

$ Comments up to 79 characters long can be used* 


Field 

Format 

Name 

Columns 

Definition 

1 

A1 

$ 

1 

Mnemonic name 

2 

A79 

COMMENT 

2-80 

Comment . 


Notes: 


1. Comments can be inserted anywhere in the input data file, except in 

between a line and its continuation line. For example, if a line has a + 
(continuation record indicator) in field 10, then the proceeding line 
cannot be commented out. 
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2.2 MMLT Data Deck 


+ 


Description: 


The + symbol indicates a continuation of the current record onto the next line. This symbol is 
placed at the end of the current line (placed in Field 10 or in column 73-80) and at the beginning of 
the subsequent line (placed in Field 1 or in column 1-8) to activate the continuation of the current 
line. 


Format and Example: 

+ 

! 2 3 4 5 6 7 8 9 

+ 

1 —2 3 4 5 6 7 8 9 

TIME 1 0.0 10.0 20.0 30.0 40.0+ 

1 2 3 4 5 6 7 8 9 

+ 50.0 60.0 


Current Line: 


Field Format 

Name 

Columns 

Definition 

10 A1 

+ 

73 

Mnemonic name 

Subsequent Line: 

Field Format 

Name 

Columns 

Definition 

1 A1 

+ 

1 

Mnemonic name 


Notes: 

1. A blank in field ten will indicate that no continuation record is to 
follow. 
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2.2 MMLT Data Deck 


2.2.1. Fabrication Process and Thermal-Mechanical Loading Data 

The following data records describe the input needed to simulate the fabrication process and 
subsequent thermal-mechanical loading of the laminate composite. The simulation of the load 
history to determine the thermal-mechanical response is represented by linear segments of varying 
temperature and mechanical loadings. The PROCESS, TIME, and TEMPR are needed in the input 
deck, the remaining mechanical loading records are optional. A list of commands in this section are 
as follows: 


Record Title 

Description 

Page 1 

Mandatory: 

PROCESS 

Defines analysis load history with linear segments. 

2 - 2 - 5 

TIME 

Specifies the time during the analysis for each linear segment. 


TEMPR 

Describes the temperature history of an individual ply. 

ESI 


Optional: 


PRESR 

Global lower and upper surface pressure history. 

2.2-8 

NX | 

Global in-plane force applied in the x-direction. 

2.2-9 

NY 

Global in-plane force applied in the y-direction. 

2.2-10 

NXY 

Global in-plane shear force. 

2.2-11 

MX 

Global moment resultant applied in the x-direction. 

2.2-12 

MY 

Global moment resultant applied in the y-direction. 

2.2-13 

MXY 

Global twisting moment resultant. 

2.2-14 


Metal Matrix Laminate Tailoring 
NASA Lewis Research Center 


Part II - MML T Input Preparation 

2 . 2-4 
















2.2 MMLT Data Deck 


PROCESS 


Description: 

The PROCESS record defines the total number of discrete load (control) points for the fabrication 
process and thermal-mechanical (TM) loads. The control points define NPNT-1 load segments. 
Subincrements are also included to specify the number of steps between control points needed for 
the nonlinear analysis of the composite. 


Format and Example : 

PROCESS NPNT I, I 2 . . . *i • • i npnt-i 

1 2 3 4 5 6 7 8 9 

PROCESS 8 5 515 20 8 10 5 


Field 

Format 

Name 

Columns 

Definition 

1 

A8 

PROCESS 

1-8 

Mnemonic name. 

2 

18 

NPNT 

9-16 

Number of load history points 
needed to define the fabrication 
process and subsequent TM load. 

3 

18 

I 

17-72 

Quantity of subincrements into which 
a segment is divided for the 


analysis. 


Notes: 

1. A maximum of 10 load points (NPNT<10) can be defined. 

2. A continuation record must be used when NPNT>7. 

3. The maximum number of subincrements is 9 or NPNT-1. 

4. The subincrement I is sequentially ordered in the PROCESS record so that 
Ij corresponds to times t ; to t i+l/ in the TIME record. 

5. In the above example NPNT=8, therefore, a total of 7 subincrements are 
needed, I, to I 7 should contain quantities. In the second segment for 
I 2 , corresponding to times t 2 to t 3 , 5 steps are assigned to perform the 
linear step analysis. 

6. The values of I, to I NPNr ., can be any integer value, but as a general 
rule in choosing the number of subincrements (I;), select I ; such that 
each subincrement represents a change in temperature of about 35°F. For 
example (refer to the TEMPR record), T i+ , = 350. 0°F and T s = 70°F then I; = 
(350-70) /35=8 subincrements. 
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2.2 MMLT Data Deck 


TIME 


Description: 

The TIME record describes time increments at each point in time for all load segments during the 
fabrication process and/or thermal-mechanical loading history. 


Format and Example: 


TIME J t j t 2 * • • 


CN • ( 
1 

| • 

1 

1 

1 

1 

tH 

1 

1 

1 

1 

1 

1 

1 + 

3 

f NPNT 
A 

i 

i 

i 

i 

i 

i 

i 

00 

1 

1 

1 

1 

1 

1 

1 

r- 

l 

l 

l 

1 

l 

l 

l 

VO 

1 

1 

1 

1 

1 

1 

1 

in 

l 

1 

l 

l 

TIME 


1 

A — 

J W ' w 

0.0 1800.0 3600.0 5400.0 5700. 

+ 6200.0 

6700.0 

7200.0 


1st Record: 




Field Format 

Name 

Columns 

Definition 

1 A8 

TIME 

1-8 

Mnemonic name. 

2-3 


9-24 

Intentionally left blank. 

4 18 

J 

25-32 

Numerical label. 

5-9 F8 

t 

33-72 

Time increments. 

10 A8 

+ 

73-80 

Continuation designation. 

2nd Record: 




Field Format 

Name 

Columns 

Definition 

1 A8 

+ 

1-8 

Continuation designation. 

2-6 F8 

t 

9-48 

Time increments. 

Notes : 





1 . The number of time steps are equivalent to NPNT specified by the PROCESS 
record, 

2. t, to t 2 are the start and end times for the first subincrement 
located in the PROCESS record at I,. 

3. The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. A + is placed in the first field of 
the second line and, beginning with the second field, the remaining time 
steps are recorded using F8 format. 

4. Time units are entered as seconds. 
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2.2 MMLT Data Deck 


TEMPR 


Description: 

The TEMPR record defines the temperature history of the individual ply corresponding to the 
number of load segments defined by the PROCESS record. 


Format and Example: 


TEMPR J Tj T 2 . . . 



-1 

2 

3 

4 — 

5- 

6- 

7 — 

8 — 

g 

+ 

-1 



3 

^NPKT 

5- 

6- 

7— 

--8 — 

9 

TEMPR 

-1 

2 

3 

2 

1500.0 
5- 

1200.0 
6- 

1075.0 
7— 

650.0 

70.0+ 

9 


+ 350.0 700.0 700.0 


1st Record: 


Field 

Format 

Name 

Columns 

Definition 

1 

A8 

TEMPR 

1-8 

Mnemonic name. 

2-3 



9-24 

Intentionally left blank. 

4 

18 

J 

25-32 

Numerical label. 

5-9 

F8 

T 

33-72 

Temperature history. 

10 

A8 

+ 

73-80 

Continuation designation. 

2nd Record: 




Field 

Format 

Name 

Columns 

Definition 

1 

A8 

+ 

1-8 

Continuation designation. 

2-6 

F8 

T 

9-48 

Temperature history. 


Notes : 

1. A temperature history must be provided for each ply in the laminate, 
therefore, a TEMPR record exists for each ply. 

2. The number of temperature loads are equivalent to NPNT specified by the 
PROCESS record. 

3. The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. A + is recorded in field 1 of the 
second line and beginning with field 2 the remaining temperature loads 
are recorded. 

4. Units for temperature are degrees Fahrenheit. 

5. Temperature T* to T i+1 is the beginning and ending temperature 
corresponding to subincrement Ij. 
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2.2 MMLT Data Deck 


PRESR 


Description: 

The PRESR command defines the pressure load history for the fabrication process of the metal 
matrix composite. Pressure is applied in the global z direction on the upper and lower surfaces. 


Format and Example: 

PRESR 

! 2 3 

■f • • * 

1 2 3 

PRESR 

1 2 3 

+ 0.0 0.0 


J 

4 


Pnpnt 
-4 — 
3 

-4-- 


0.0 


Pj P2 • • • 

— 5 6 7 8 

— 5 6 7 8 

0.0 - 1000.0 - 1000.0 - 1000.0 

— 5 6 7 8 


Pi+ 

9 


— 9 
0 . 0 + 
— 9 


1st Record: 


Field 

Format 

Name 

Columns 

1 

A8 

PRESR 

1-8 

2-3 



9-24 

4 

18 

J 

25-32 

5-9 

F8 

P 

33-72 

10 

A8 

+ 

73-80 

2nd Record: 



Field 

Format 

Name 

Columns 

1 

A8 

+ 

1-8 

2-6 

F8 

P 

9-48 


Definition 
Mnemonic name. 
Intentionally left blank. 
Numerical label* 

Pressure history. 
Continuation designation. 

Definition 

Continuation designation. 
Pressure history. 


Notes: 


1. The number of pressure loads are equivalent to NPNT specified by the 
PROCESS record. 

2. The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. A + is placed in the first field of 
the second line and beginning with the second field the remaining 
pressure loads are recorded. 

3. Units for pressure are pounds per square inch (psi) . 

4. Pressure Pj to P i+t is the beginning and ending pressure corresponding to 
subincrement I,. 

5. This command is optional in the input deck. 
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2.2 MMLT Data Deck 


NX 


Description: 

A force applied to the laminate composite in the global x direction. Forces are applied at each load 
point with the default equaling zero. 


Format and Example: 

NX J N xt N*2 . . . 


+ 

i 

04 • C 

1 

| • 

1 

| ♦ 

1 

1 

1 

iH f 

1 


N xNPNT 



NX 

— x z 


3 4 

4 

0.0 0.0 0.0 0.0 

0. 

+ 

“ 1 2 

0.0 

0 . 

0 500.0 



1st Record: 





Field 

Format 

Name 

Columns 

Definition 


1 

A8 

NX 

1-8 

Mnemonic name. 


2-3 



9-24 

Intentionally left blank. 


4 

18 

J 

25-32 

Numerical label. 


5-9 

F8 

N x 

33-72 

Applied resultant forces. 


2nd Record: 





Field 

Format 

Name 

Columns 

Definition 


1 

A8 

+ 

1-8 

Continuation designation. 


2-6 

F8 

N, 

9-48 

Applied resultant forces. 



Notes : 

1* The number of force loads are equivalent to NPNT specified by the 
PROCESS record. 

2. The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. Refer to the continuation record (+) 
for further details. 

3. Forces N^ to N^+i are the beginning and ending force corresponding to 
subincrement and vary linearly from tj to t i+l . 

4. The force has units of force per unit depth. 

5. A blank field constitutes a zero entry. 

6. This command is optional in the input deck. 
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2.2 MMLT Data Deck 


NY 


Description: 

A force applied to the laminate composite in the global y direction. Forces are applied at each load 
point with the default equaling zero. 


Format and Example: 
NY 


N v 


■V 


+ 

-1 2- 

3 

4 

NyNPNT 

A 

1 Ul 

> <J> 

J '-J 

NY 

— 1 2- 


5 

A _ 

W w * 

0.0 0.0 0.0 < 

+ 

— 1 2- 

0.0 

0.0 

500.0 


1st Record: 




Field 

Format 

Name 

Columns 

Definition 

1 

A8 

NY 

1-8 

Mnemonic name. 

2-3 



9-24 

Intentionally left blank. 

4 

18 

J 

25-32 

Numerical label. 

5-9 

F8 

N y 

33-72 

Applied resultant forces. 

2nd Record: 




Field 

Format 

Name 

Columns 

Definition 

1 

A8 

+ 

1-8 

Continuation designation. 

2-6 

F8 

Ny 

9-48 

Applied resultant forces. 


- 8 - 

- 8 - 
0.0 
8 - 


Nyi+ 


9 

0 . 0 + 

9 


Notes: 


1. The number of force loads are equivalent to NPNT specified by the 

PROCESS record. . . . . 

2 The continuation record is needed when NPNT>5. This is shown in the 

above Format and Example section. Refer to the continuation record (+) 
for further details. 

3. Forces N,* to N yi+ , are the beginning and ending force corresponding to 
subincrement Ij and vary linearly from t; to t i+) . 

4. The force has units of force per unit depth. 

5. A blank field constitutes a zero entry. 

6. This command is optional in the input deck. 
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2.2 MMLT Data Deck 


NXY 


Description: 

A shear force applied to the laminate composite in the global xy plane. The shear forces are applied 
at each load point with the default equaling zero. 


Format and Example: 


NXY 




J 

Nxy, 


• . 


Nxyi + 


— 1— 

2 

3 — 

4 

5 

6 

7 

8 

9 

+ 


• • • 


NjyNPNT 




g 



— 1 — 

2 

3 — 

4 

5 

— — 6— — — - 

7 

— 9 

NXY 




6 

0.0 

0.0 

0.0 

0 . 0 

0.0+ 


— 1 — 

2 

3 — 

4 

-———5 

6- 

7 

8— 

9 

+ 


0.0 

0 . 0 

500.0 







1st Record: 


Field 

Format 

Name 

Columns 

Definition 

1 

A8 

NXY 

1-8 

Mnemonic name. 

2-3 



CN 

1 

crv 

Intentionally left blank. 

4 

18 

J 

25-32 

Numerical label. 

5-9 

F8 


33-72 

Applied shear forces. 


2nd Record: 


Field 

Format 

Name 

Columns 

Definition 

1 

A8 

+ 

1-8 

Continuation designation 

2-6 

F8 

N *y 

9-48 

Applied shear forces. 


Notes: 

1. The number of shear forces are equivalent to NPNT specified by the 
PROCESS record. 

2. The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. Refer to the continuation record ( + ) 
for further details. 

3. Forces N^ to N^, are the beginning and ending force corresponding to 
subincrement I* and vary linearly from t 5 to t i+1 . 

4. The force has units of force per unit depth. 

5. A blank field constitutes a zero entry. 

6. This command is optional in the input deck. 
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2.2 MMLT Data Deck 


MX 


Description: 

A moment applied to the laminate composite in the global x direction. Moments are applied at each 
load point with the default equaling zero. 


Format and Example: 


NX 




J 

5 — - 

6 * . 

* • _ 

8 


+ 

- — 1 — 



MxNPNT 

5— 



8 

9 

NX 

1 — 

1 — 



7 

0.0 

5 

0.0 

0.0 

0.0 

8 

0 . 0 + 
9 


0.0 0.0 10.0 


1st Record: 




Field 

Format 

Name 

Columns 

Definition 

1 

A 8 

MX 

1-8 

Mnemonic name. 

2-3 



9-24 

Intentionally left blank. 

4 

18 

J 

25-32 

Numerical label. 

5 

F 8 


33-72 

Applied resultant moment. 

2nd Record: 




Field 

Format 

Name 

Columns 

Definition 

1 

A 8 

+ 

1-8 

Continuation designation. 

2-6 

F 8 


9-48 

Applied resultant moment. 


Notes: 

1. The number of moment loads are equivalent to NPNT specified by the 

PROCESS record. ... . . 

2. The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. Refer to the continuation record (+) 
for further details. 

3. Moments M*; to M^, are the beginning and ending moment corresponding to 
subincrement I s . 

4. The moment has units of lbs-in/in. 

5. A blank entry constitutes as a zero entry. 

6 . This command is optional in the input deck. 
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MY 


Description: 

A moment applied to the laminate composite in the global y direction. Moments are applied at each 
load point with the default equaling zero. 


Format and Example: 

NY J My! My 2 . • • 


+ 

— 1— 

* • • 


My NPNT 

5 

6 



NY 

1 — 



■^4 

8 

5 ■ 

0.0 

0.0 

0.0 

0.0 

+ 

1 — 

0.0 

— — 2 

0.0 

5.0 






1st Record: 




Field 

Format 

Name 

Columns 

Definition 

1 

A8 

MY 

1-8 

Mnemonic name. 

2-3 



9-24 

Intentionally left blank. 

4 

18 

J 

25-32 

Numerical label. 

5-9 

F8 


33-72 

Applied resultant moments 

2nd Record: 




Field 

Format 

Name 

Columns 

Definition 

1 

A8 

+ 

1-8 

Continuation designation. 

2-6 

F8 

My 

9-48 

Applied resultant moment. 


-9 


— 9 
0 . 0 + 
— 9 


Notes: 

1. The number of moment loads are equivalent to NPNT specified by the 
PROCESS record. 

2. The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. Refer to the continuation record (+) 
for further details. 

3. Moments to M yi+l are the beginning and ending moment corresponding to 
subincrement I; and vary linearly from t t to t i+1 . 

4. The moment has units of lbs-in/in. 

5. A blank entry constitutes as a zero entry. 

6. This command is optional in the input deck. 
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MXY 


Description: 

A moment applied to the laminate composite in the global xy plane. The moment are applied at 
each load point with the default equaling zero. 


Format and Example: 
NXY 


+ 
NXY 


- 1 - 

- 1 - 


-3- 

-3- 


J 

—4 

^xyNPNT 

— 4 


H*, 

5- 

5- 

0.0 


6 — 

6 — 

0.0 


7- 

0.0 


! 2 


•3 4- 


+ 0.0 

0. 

0 3.0 


1st Record: 

Field Format 

Name 

Columns 

Definition 

1 A8 

MXY 

1-8 

Mnemonic name. 

2-3 


9-24 

Intentionally left blank. 

4 18 

J 

25-32 

Numerical label. 

5-9 F8 


33-72 

Applied resultant moment. 

2nd Record: 

Field Format 

Name 

Columns 

Definition 

1 A8 

+ 

1-8 

Cont inuat ion des ignat ion . 

2-6 F8 


9-48 

Applied resultant moment. 

Notes: 

1. The 

number 

of moment 

loads are equivalent to NPNT 


— 8 - 

— 8 - 
0.0 
8 - 


9 

9 

0 . 0 + 

9 


4. 

5, 
6 « 


PROCESS record. 

The continuation record is needed when NPNT>5. This is shown in the 
above Format and Example section. Refer to the continuation record (+) 
for further details. 

Moments to are the beginning and ending moment corresponding to 

subincrement I* and vary linearly from tj to t i+1 . 

The moment has units of lbs-in/in. 

A blank entry constitutes as a zero entry. 

This command is optional in the input deck. 
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2.2.2 Material and Laminate Identification 

In this section the laminate make-up, constituent materials and ply characteristics are identified for 
the desired laminate composite. The following list identifies the records in this section. 



Mandatory: 


LAY-UP 

Defines the lamina order for a composite material. 

2.2-17 

MATERIAL 

Defines the constituent and ply characteristics. 

2.2-18 
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PRIORITY 


Description: 

The PRIORITY record defines the hierarchy used to assign the ply thicknesses for the laminate. If 
the PRIORITY record is omitted the LAY-UP record input will define the layer hierarchy. 


Format and Example: 

PRIORITY plyj ply 2 • * • 

1 2 3 4* 

PRIORITY SBL1 SBL2 SBL3 


plyj . . • piy* 

5 6 7 8 9 


Field Format Name 


Columns Definition 


1 

A8 

PRIORITY 

1-8 

2-9 

A8 

ply 

9-72 


Mnemonic name. 

MATERIAL names listed in descending 
order of importance. 


Notes : 


1 . 

2 . 


3. 


4. 

5. 


A maximum of 10 plies can be used. 

This card should only be used when the thickness is used as a design 

variable. . . 

In the above example, the ply# named SBL1# would be given highest 
priority, due to the layer hierarchy, followed by each successive layer 
in the list. The entire thickness of layer SBL1 would be used and if 
the total thickness was not filled then the SBL2 layer would be added. 
This process would continue until the total thickness was filled. 

The ply values are the same values indicated in the MATERIAL records. 
This record is optional in the input deck. 
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LAY-UP 


Description: 

The LAY-UP record defines the ply lay-up for the specific composite listing of each individual ply 
starting at the bottom ply and moving upward through the thickness of the laminate. 


Format and Example: 


LAY-UP ply, ply 2 - • • P l Y» 

1 2 3 4 5 6 7 8 9 

LAY-UP SBL1 SBL2 SBL3 


Field Format Name 

1 A8 LAY-UP 

2-9 A8 ply 


Columns Definition 

1-8 Mnemonic name. 

9-72 MATERIAL names listed from bottom 

surface to top surface. 


Notes: 

1. A maximum of 10 plies can be used. 

2. All plies should be listed in their proper order, even for the case of a 
symmetric lay-up. 

3. The ply names are used to identify each separate ply and the same names 
are utilized in the MATERIAL records. 

4. In the example, SBL1, SBL2 , and SBL3 represent the individual plies 1, 

2, and 3, respectively. 

5. Refer to Figure 1.5-3 for ply ordering sequence. 
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MATERIAL 


Description: 

The MATERIAL record defines the ply properties of the individual plies that make-up the laminate, 
for example, the fiber volume ratio (FVR), interphase thickness, ply thickness, and the ply 
orientation. 


Format and Example: 


MATERIAL 

ply 

CODE 

FVR 

A 


MATERIAL SBL1 PIOOCOPR 0.40 0.12 


Field 

Format 

Name 

Columns 

1 

A8 

MATERIAL 

1-8 

2 

A8 

ply 

9-16 

3 

A8 

CODE 

17-24 

4 

F8 

FVR 

25-32 

5 

F8 

t„ 

33-40 

6 

F8 


41-48 

7 

F8 

ANGLE 

49-56 


t, ANGLE 

6 7 8 9 

0.01 0.0 

Definition 
Mnemonic name. 

User supplied name for ply. 

Fiber and matrix code names. 

Fiber volume fraction. 

Interphase thickness. 

Ply thickness. 

Angle orientation (degrees) relative 
to the global axis. 


Notes : 

1. A separate MATERIAL record must be defined for each ply in the 
composite. 

2. A maximum of 10 different plies can be specified. 

3. The ply name (e.g. SBL1) should match a name given in the PRIORITY and 
LAY-UP records and should be unique to all other plies. 

4. The CODE field specifies the four letter code names of the fiber and 
matrix needed to identify the corresponding entries in the constituent 
databank [Reference 3]. A description on the setup for the databank is 
located in Appendix B. The four letter code name for each the fiber and 
matrix is specified in the databank and matched in the CODE field. 

5. In the above example, the CODE field specifies PIOOCOPR for a 
graphite /copper composite, therefore, in the databank file P100 would 
specify the constituent properties for the fiber and the COPR specifies 
the copper matrix properties. 

6. The user has the ability to setup their own code names within the 
databank. 
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7. The fiber volume fraction is restricted to the range: 

0.05 < FVR < 0.78 

8. If the interphase thickness is: 

(i) t d > 0.0 An interphase exists within the model and 
extends into the matrix. 

( ii ) t d = 0.0 No interphase exists within the model. 

(iii) t d < 0.0 An interphase exists within the model and grows 

into the fiber. 

The interphase is entered as a percentage of the fiber diameter. In the 
above example the interphase has a 0.12 or 12% of the fiber diameter 
thickness. 

9. Perfect bonding between the fiber and matrix is assumed if no interphase 
is considered (i.e. if t d =0.0). 

10. For the angle orientation, 0, refer to Figure 1.5-4. In the example the 
composite is an unidirectional ply, as a result, 0 is equal to 0.0. 
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2.2.3 Design Variables and Constraints 

The design variables and constraints are described by the following commands. 


1 Record Title 

Description 

Page 

II VARIABLE 

Defines an individual design variable. 

2.2-20 

DEPEND 

Associates a dependent variable to a design variable. 

2.2-22 

CONSTRNT 

Defines the specific constraints for the optimizer. 

2.2-23 


Metal Matrix Laminate Tailoring 
NASA Lewis Research Center 


Part II - MML T Input Preparation 

2.2-20 










2.2 MMLT Data Deck 


VARIABLE 


Description: 

The VARIABLE record defines the information for an individual design variable (parameter) used 
for the tailoring procedure. 


Format and Example: 

VARIABLE name J VPOS Z L Z u Z_START 

VARIABLE TEMPI 1 11 70.0 1500.0 1350.0 


Field 

Format 

Name 

Columns 

Definition 

1 

A8 

VARIABLE 

1-8 

Mnemonic name. 

2 

A8 

name 

9-16 

User supplied name for variable. 

3 

18 

J 

17-24 

Numerical identification label. 

4 

18 

VPOS 

25-32 

Position of variable within vector 
of available variables. 

5 



33-40 

Intentionally left blank. 

6 

F8 

Z L 

41-48 

Design variable lower bound. 

7 

F8 

z u 

49-56 

Design variable upper bound. 

8-9 

F16 

Z_START 

65-72 

Initial value of design variable. 


Notes : 

1. Available design variables and their locations are identified in Table 

2 . 2 - 1 . 

2. Each design variable has an individual VARIABLE record. For example , if 
there are 5 design variables then there must be 5 VARIABLE cards. 

3. Initial values of the design variables will overwrite all similar data 
given by previous input commands. For example, if the consolidation 
pressures are given in the PRESR command and are also design variables, 
the initial design variables will overwrite the corresponding pressures 
given in the PRESR record. 

4. In the example, TEMPI represents the temperature of ply 1 for the first 
load step which is defined by the VPOS=ll. 

5. VPOS points to the location of the design variable in the variable 
storage vector. Further details are located in Table 2.2-1 on the 
following page. 

6. The numerical label (J) should be unique for each VARIABLE record used 
within the input deck and defined in a numerical order. 
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Table 2.2-1: Mapping of Available Design Variables 


This table lists the available design variables, their dimensions, and locations 
Si£hin the storage array. The user is referred to the VARIABLE record field 4 
(VPOS) and the DEPEND record fields 3 and 4 (VPOSl and VPOS2 ) , respective y. 

— 

Description 

Dimension 

Position 

Time 

NPRS 

1-10 

Temperature 

NPRS*NLR 

11-110 

Pressure 

NPRS 

111-120 

Fiber Orientation 

NLR 

121-130 

Fiber Volume Ratio 

NLR 

131-140 

Interphase Thickness 

NLR 

141-150 

Interphase Modulus 

NLR 

151-160 

Interphase CTE 

NLR 

161-170 

Interphase Strength 

NLR 

171-180 

Zeta 

1 

181 

Force , N x 

NPRS 

191-200 

Force, N y 

NPRS 

201-210 

Force , Njjy 

NPRS 

211-220 

Moment , M* 

NPRS 

221-230 

Moment , My 

NPRS 

231-240 

Moment , M^ 

NPRS 

241-250 


Notes : 


1. NPRS is the maximum number of load history segments. Refer to the 
PROCESS record: 

(i) NPRS = 10 

(ii) NPNT < NPRS 

2. NLR is the maximum number of plies the analysis can model. 

(i) NLR =10 

3. The VARIABLE and DEPEND records define the storage position 
VPOSl, and VPOS2 ) of the available design parameters. 

4. Zeta is the upper bound for the min/max objective function, 

5 For T exLiple? r if the ply orientations (for a laminate with 3 plies) are 

specified as the design variables, then VPOS for ply 1 (bottom ply) 
equal to 121, likewise, for ply 2 VPOS = 122 and ply 3 (top ply) VPOS = 
123. 


(VPOS, 
see the 
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DEPEND 


Description: 

The DEPEND record links additional design parameters which are dependent upon a specified 
VARIABLE. The DEPEND record allows the user to alter input values without defining another 
independent design variable. A multiplier is permitted. 


Format and Example: 


DEPEND 

name 

VPOSl 

o _ 

VPOS2 RATIO 

A C 


DEPEND 
Field : 

1 z — 

Format 

21 

Name 

11 1.0 
Columns 

Definition 

1 

A8 

DEPEND 

1-8 

Mnemonic name . 

2 

A8 

name 

9-16 

User supplied name. 

3 

18 

VPOSl 

17-24 

Position of variable within storage 

4 

18 

VPOS 2 

25-32 

vector of available variables. 

Location (VPOS) of the VARIABLE design 

5 

F8 

RATIO 

33-40 

value that the additional variable is 
dependent upon. 

Value assigned as multiplier. 


Notes: 

1. Positions (VPOS1 and VPOS2 ) for variables are located in Table 2.2-1. 

2. VPOS2 is equivalent to VPOS in the VARIABLE record that the DEPEND 
record is linking to. 

3. In the above example, refer to Table 2.2-1, the DEPEND records indicates 
that the temperature for the second ply at the first load point, defined 
by VPOSl=21, will be changed at a 1:1 ratio (RATIO=1.0) to the design 
variable that is defined by VPOS2=ll. A VARIABLE record must exist 
defining the temperature at storage location (VPOS) 11 which indicates 
that the temperature in the first ply at the first load point is a 
design parameter. 
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CONSTRNT 


Description: 

The CONSTRNT record defines the constraints for the optimization program. All information 
needed to compare the calculated value to the allowable value is given by this record along with the 
specified equation form (refer to Table 2.2-2). The generated terms are stored in an array and 
referenced by array location (refer to Table 2.2-3). 


Format 

and Example 

S 




CONSTRNT NAME 

FACTOR 

I FORM ITYPE 

A C 

ITERM1 ITERM2 

&CNST 

-1 2“ 

VALl 

VAL2 

A 

_c_ 


CONSTRNT STRESS 

1.0 

1 

A _ . 

0 

151 158 

&CNST 

-1 2 — 

0.95 





1st Record: 





Field 

Format 

Name 

Columns 


Definition 

1 

A8 

CONSTRNT 

1-8 


Mnemonic name. 

2 

A8 

NAME 

9-16 


User supplied constraint name. 

3 

F8 

FACTOR 

17-24 


Multiplication factor applied to 
TERMS for analysis calibration. 

4 

18 

I FORM 

25-32 


Equation form of constraint. 

5 

18 

ITYPE 

33-40 


0 ■ nonlinear inequality constraint 
2 = linear inequality constraint 
-1 = nonlinear equality constraint 

6 

thru 

7 

18 

ITERM1 

thru 

ITERM2 

41-56 


Location of constraint terms 
in storage array, TERMS. 

2nd Record: 





Field 

Format 

Name 

Columns 


Definition 

1 

A8 

&CNST 

1-8 


Mnemonic name. 

2 

2F8 

VALl 

thru 

VAL2 

9-24 


Value needed to calculate constraint 
equations listed in Table 2.2—2. 
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Notes: 

1. The CONSTRNT record must have the 2nd record. 

2. Refer to Table 2.2-2 to set up the proper constraint equations. 

3. Available constraints and their locations are identified in Table 2.2-3 
to properly set up constraint equations. 

4. TERMS refers to the actual values stored in the constraint arrays and 
located by the ITERM pointers. 

5. ITERM1 and ITERM2 define the initial and final locations in TERMS array. 
In the above example, eight constraints are defined. In case of single 
constraint definitions ITERM2 may be omitted. 


Table 2.2-2: Constraint Equations 


The following table is a list of available constraint equations that can be 
used for the optimizer in the MMLT program. The ID Number listed before each 
equation is the input on the CONSTRNT record for field 4, IFORM. Each 
equation compares one or more of the VAL1 through VAL2 input values to the 
constraint terms, located by ITERM1 through ITERM2 on the CONSTRNT record 
(refer to Table 2.2-3). Equations are interpreted as follows: 

G < 0 inactive constraint 
G = 0 active constraint 
G > 0 violated constraint 

ID 

Number 

Form of Equation 

Description 

1 

G = (TERM - VAL1) / VAL1 

Upper bound constraint. 

2 

G = (VAL1 - TERM) / VAL1 

Lower bound constraint. 

3 

G = VAL2 - (TERM - VAL1) / VAL1 

Constraint from above with margin. 

4 

G = (TERM - VAL1 ) / VAL1 - VAL2 

Constraint from below with margin. 
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This table defines the storage locations of the. specified 
constraints calculated during the optimization. The positrons 
are given under the appropriate heading and are input on the 
CONSTRNT record fields 6 and 7, ITERMS. 


Description 

(TERMS) 


Dimension 


Zeta,a rallA 


Zeta , a 


Constraints ( Imposed at each control point 

NPRS*NLR 


NPRS*NLR 


NPRS*NLR 


NPRS*NLR 


Position 

(ITERMS) 



51-150 


151-250 


251-350 


351-450 


451-550 


551-650 


651-750 


751-850 
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u> to 
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Description 

Dimension 

Position 

(TERMS) 


(I TERMS) 

St22 

NLR 

1871-1880 

Sc22 

NLR 

1881-1890 

s, 2 

NLR 

1891-1900 

S,3 

NLR 

1901-1910 

S23 

NLR 

1911-1920 

1 

Interphase Thickness 

NLR 

1921-1930 

Extensional Stiffness 

3 

1931-1933 

|| Flexural Stiffness 

3 

1934-1936 

Temperature 

NPRS 

1937-1946 

Pressure 

NPRS 

1947-1956 


Notes: 


1 . 


4. 


5. 


NPRS is the maximum number of load history steps (10) and corresponds to 
the NPNT of the PROCESS record. 

NLR is the maximum number of plies (10) in the composite. 

For example, if NPNT=8 and the composite has 3 plies and the fiber 
stress constraint, o m , is used. A CONSTRNT record is needed for 
each ply and is represented by the following positions in the storage 
array TERMS : 

a. ITERM1 = 751 and ITERM2 = 758, for ply 1 

b. ITERM1 = 761 and ITERM2 = 768, for ply 2 

c. ITERM1 = 771 and ITERM2 » 778, for ply 3 

A constraint should only be used when it makes physical sense, for 

example, if no interphase is used, then stress failure constraint in the 
B region of the composite model should not be included in the activated 
constraints. 

The stress failure constraints are represented by the following 
equation: 

S c < a < S T 

which is applied throughout the loading history. 
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2.2.4 Optimization Parameters and Objective Function 

This section describes the input records required to set up the optimization problem and the 
objective function solution. The following records are explained in this section. 


Record Title 

Description 

Page 

OPTIMIZE 

Sets up the optimization strategy and procedure. 

2.2-28 

OBJTYPE 

Objective function equation type. 

2.2-29 

SENSE 

Sensitivity values for the objective function. 

2.2-30 

*END_INP 

Signals the end of input deck. 

2.2-32 


Metal Matrix Laminate Tailoring 
NASA Lewis Research Center 


Part II - MMLT Input Preparation 

2.2-28 








2.2 MMLT Data Deck 


OPTIMIZE 


Description: 

The OPTIMIZE record defines the optimization solution procedure, the optimization strategy, and 
the search procedure. 


Format and Example: 


OPTIMIZE 

ISTRAT 

IOPT 

ISERCH 

IOUT 

ISCAL 

IDM 

MXITER 

!- 

2 — 

3- 

4 — 

5— 

—————— 6— — 

7- 

8— 

OPTIMIZE 

0 

5 

8 

3440 

1 

40 

1000 


Field 

Format 

Name 

Columns 

Definition 

1 

A8 

OPTIMIZE 

1-8 

Mnemonic name. 

2 

18 

ISTRAT 

9-16 

Optimizer strategy. 

3 

18 

IOPT 

17-24 

Optimization procedure. 

4 

18 

ISERCH 

25-32 

Optimizer search procedure. 

5 

18 

IOUT 

33-40 

Optimizer output request. 

6 

18 

ISCAL 

41-48 

Optimizer scaling indicator. 

7 

18 

IDM 

49-56 

Desired number of completed 
optimizer design moves. 

8 

18 

MXITER 

57-64 

Maximum number of objective 
function/constraints evaluations 


Notes: 

1. The user should use the values given in the example for input. The MMLT 
code was not hardwired for these values to give the user flexibility if 
they chose to change any of these values. 

2. For further information on the terms in this record, the user is 
referred to Reference 4. 

3. One optimization step requires at least NDV+1 objective 

function/ constraint evaluations, where NDV is the total number of design 
variables defined by the VARIABLE records. 
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OBJTYPE 

Description: 

The OBJTYPE record defines the objective function equation type. 




Ob j active Function (OF) 


Description 


OF * S,*zeta 


2 OF = S,*zeta + S 2 *(T,/T 0 ) 


OF = S,* (o m]1A ) 2 + S 2 *(«W) + S 3 *(o 01 33 A ) 2 

+ S^Jo^.b) 2 + S 5 *(o,»22b) 2 + S 4 * (o^jjb) 2 
+ S 7 *(<W) 2 + S ,,*(CW)* + S„*(g m3 3c) 


OF = (S,*NJ + (S 2 *N y ) + (Sj^Njy) 

+ (S 4 *MJ + (Sj*^) + (S^) 

+ <S 7 *P) + ( S g *T ) 


OF = (S,*A„) + (S 2 *A22) + (S 3 *Atf) 
+ (S 4 *D„) + (S 5 *D h ) + (S^D*) 


minimize/maximize design 
variable zeta. ____ 


minimize/maximize design 
variable zeta with temperature 
added . 


minimize/maximize specific 
matrix stresses. 


minimize/maximize specific 
thermal-mechanical loads . 


minimize/maximize the laminate 
stiffness, extensional or 
bending* 


where Sj is the weighted factor. 
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SENSE 


Description: 

The SENSE record defines the sensitivity factors for the objective function shown in the OBJTYPE 
record. The sensitivity factors are used as weighted coefficients to determine rank of importance. 


Format and Example: 


SENSE 

S, 

s 2 

s 3 S 4 

4 5 

s s 

s 6 s 

SENSE 

— 1 2- 

1.0 

1.0 

0.0 0.0 

0.0 

0.0 0.1 

1st Record: 





Field 

Format 

Name 

Columns 

Definition 

1 

A8 

SENSE 

1-8 

Mnemonic 

name. 

2 

F8 

S, 

9-16 

Weighted 

coefficient . 

3 

F8 

s 2 

17-24 

Weighted 

coefficient . 

4 

F8 

S 3 

25-32 

Weighted 

coefficient. 

5 

F8 

S 4 

33-40 

Weighted 

coefficient. 

6 

F8 

S 5 

41-48 

Weighted 

coefficient* 

7 

F8 

S« 

49-56 

Weighted 

coefficient. 

8 

F8 

S 7 

57-64 

Weighted 

coefficient. 

9 

F8 

^8 

65-72 

Weighted 

coefficient . 

2nd Record: 





Field 

Format 

Name 

Columns 

Definition 

1 

A8 

+ 

1 

Mnemonic 

name. 

2 

F8 

Sg 

9-16 

Weighted 

coefficient . 

3 

F8 

S]0 

17-24 

Weighted 

coef f icient . 
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in 


2.2 MMLT Data Deck 


SENSE (cont.) 


Notes: 

1. If the sensitivity is equal to 

(a) a positive value the respective term will be minimized 

(b) a negative value the respective term will be maximized 
based on a positive design value. 

2. The values S, thru S, 0 correspond to the objective functions listed under 
the OBJTYP record and the definition of the ITYP value. In the above 
example, if ITYP=4, then S, and S 2 activate and equally weighted the N x 
and N y terms in the objective function. 

3. If the sensitivity factor is activated by a nonzero number and ITYP-4 in 
the OBJTYP record, a VARIABLE record must exist for the corresponding 
load, otherwise, an error will occur. 

The continuation card is optional. 

Default value is 0.0. 
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*END_INP 

Description: 

The *END_INP record signals the end of the input file. 

Format and Example: 

*END INP 

Z — i 2 3 4 5 6 7 8 9 

*END_INP 

Field Format Name Columns Definition 

1 A8 *END INP 1-8 Mnemonic name. 
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Two examples of input data files, shown in Tables 2.3-1 and 2.3-2, are presented in this section to 
aid the user in setting up an input file for the MMLT computer code. In addition, the input files 
demonstrate the flexibility and capabilities of the methodology presented in this manual. In the first 
example, a [0/90] 2S graphite (P100)/copper (Cu) MMC is concurrently tailored for the fabrication 
process parameters and interphase characteristics in order to reduce the residual matrix microstresses 
at the end of fabrication. The first example gives a more detailed description than the second 
example, with each record in the input file given an explanation. The second example involves a 
[±45]s silicon carbide (SiC)/titanium (Til5-3-3-3) with the fabrication process parameters tailored to 
maximize the in-service loading. In the following sections a description and input file are provided 
for these two examples. 


2.3.1 Case 1: [0/90J 2S PlOO/Copper 

The first example problem considers a [0/90]2s laminate made of ultra high modulus graphite 
(P100)/copper (Cu). The objective is to tailor the fabrication process parameters and interphase 
characteristics by reducing the longitudinal and transverse stresses in the matrix. Typical design 
variables include temperature and pressure histories, the interphase modulus, CTE, strength, and 
layer thickness. 

Only the cool-down phase of the fabrication process is simulated during the tailoring and the 
thermo-mechanical response of the composite. It can be assumed that residual stresses are 
negligible in the beginning phases of fabrication. The same fabrication process used for the example 
in Section 1.2 for unidirectional PlOO/Cu was used as the current process for this laminate (refer to 
Figure 1.3-1). The first few lines of the input file in Table 2.3-1 are comment records ($) used to 
describe the problem. A PROCESS record follows to indicate the load history defined by the 
control points in Figure 1.3-1. The first integer value represents the total number of control points 
(NPNT=5) and has NPNT-1 linear segments represented by the subincrements (Ii-npnt-i) between 
control points. The TIME, TEMPR, PRESR, and NY records are assigned next. These records 
simulate the cool-down phase of the fabrication process. 


To identify the composite lay-up and material characteristics, the LAY-UP and MATERIAL 
records are used. The LAY-UP record indicates the individual plies (e.g., PLY1, PLY2, and 
PLY3) in their stacking sequence by a user supplied name. A MATERIAL record contains the user 
supplied name of the ply (e.g. PLY1), composite system (PIOOCOPR), the FVR (0.40), interphase 
layer thickness (12% of the fiber diameter), ply thickness (e. g. 0.01), and angle of ply (e.g. 00.0). 
For each ply a name must be given in the LAY-UP and MATERIAL records that is the same in 

both records. 

A VARIABLE record identifies the 1 1 design variables that are tailored. These are the 
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2.3 Sample Input Data Files 


temperature history (3 design variables) where the consolidation and room temperatures are not used 
in the tailoring, the pressure history (4 design variables) where the final pressure is always zero, and 
the interphase properties (4 variables: modulus, CTE, strength, and layer thickness). The current 
process parameters are shown in Figure 1.3-1. Initial interphase properties were assumed to be the 
same as the matrix properties. For example, the VARIABLE record can be defined by a user 
supplied name (e.g., TEMPI 1), a numerical label (e.g., 1), the position of the variable in the 
storage array (12, refer to Table 2.2-1), the lower bound, the upper bound, and the initial or current 
value. 


The DEPEND records are used in this example so that consistency among the three plies and 
loading is maintained. In effect, no new design variables are introduced in tailoring a laminate due 
to its extra plies. Only the temperature history for ply 1 is used as design variables but the 
temperature histories for plies 2 and 3 are scaled to equivalent values as ply 1 . For the pressures 
applied to the upper and lower surfaces at each load point, the resultant force (N y ) applied in the y- 
direction is scaled to the pressure times laminate thickness. This keeps isostatic pressure in the y 
and z planes. The interphase characteristics are scaled in the same manner as the temperature 
history, i.e., the interphase characteristics for ply 1 are set as the design variables, where plies 2 
and 3 are scaled to equivalent interphase values that are in ply 1 . A DEPEND record needs a user 
supplied name (e.g. TEMP22), the position in the storage array (e.g. 22, refer to Table 2.2-1), the 
position of the design variable that it is dependent upon (e.g. 12), and a multiplication factor (e.g. 
1 . 0 ). 


Stress constraints are imposed at the control points and defined in the PROCESS record to 
insure the integrity of the material at all phases of the tailoring. In this manner, five constraints 
were introduced for each microstress (12) in each ply (3), therefore, a total of 180 stress constraints 
exist. The CONSTNT and &CNST records identifies the constraint name (e.g. SfllCl), 
multiplication factor (1.0), equation form (1, refer to Table 2.2-2), type of constraint (0 indicating a 
nonlinear constraint), location of constraint in storage array (refer to Table 2.2-3), and for the 
second record a value that is used to evaluate the constraint. Finally, three additional constraints on 
the stiffness of the laminate are included to ensure that the integrity of the tailored composite is 
within 85 % of the current composite. Therefore, no substantial degradation is possible during 
tailoring. 

An OPTIMIZE record is used to activate the tailoring procedure by defining the optimizer 
strategy (0), optimization procedure (5), search procedure (8), output request (3440), scaling 
indicator (1), maximum number of design moves, and the maximum number of evaluations (200). 
For further details consult Reference 4. To identify the objective function, the OBJTYPE record is 
used. In this case, the maximum matrix microstresses (a mUA and < r my >* ) are reduced at the end of 
fabrication, therefore, 3 identifies this objective function. The SENSE record is used to indicate the 
active matrix microstresses in the objective function. Finally, the *END_INP record is used to 
indicate the end of input. 
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Table 2.3-1: Input for Case 1 (PlOO/Copper [0/90]s) 


$= 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 


Material 

Tailor 


Objective 


no. of load history 


PROCESS 5 


CASE 1 

MMLT Input Data File 

: - [0/90)25 Graphite ( P100 ) /Copper (COPR) 

: - Fabrication Parameters 

* Temperature and Pressure History 
- Interphase Characteristics 

* Modulus, CTE, Strength, Thickness 

: - Minimize the Residual Matrix Microstresses at 

the End of Fabrication in the longitudinal (11) 
and transverse (22) directions 


Fabrication Process 

points (NPNT), no. of subincrements (I(NPNT-l)) 
10 40 10 5 


$ 

$ admissible parameters: 
$ time (secs) 

TIME 

$ temperature history 


TEMPR 

TEMPR 

TEMPR 


TIME, TEMPR, PRESR, NX, NY, NXY, MX, MY, MXY 

1 0.0 1801.0 3601.0 5401.0 5701.0 

(°F) of cool-down phase, a TEMPR record for each ply 

2 1500.0 1500.0 1075.0 650.0 70.0 

3 1500.0 1500.0 1075.0 650.0 70.0 

4 1500.0 1500.0 1075.0 650.0 70.0 


$ consolidation pressure history 

PRESR 5 


(psi) 

0.0 - 1000.0 - 1000.0 - 1000.0 


0.0 


$ 

NY 


resultant force (lb/in), equivalent to P; * t, 

6 0.0 -40.0 -40.0 


-40.0 0.0 


$ Material and Composite Identification 

$ assign ply lay-up 


LAY-UP 

PLY1 

PLY2 

PLY 3 




$ 

$ 


fib/mat 

FVR 

t d 

ti 

ang 

MATERIAL 

PLY1 

P100COPR 

0.40 

0.12 

0.01 

00.0 

MATERIAL 

PLY 2 

P10QCOPR 

0.40 

0.12 

0.02 

90.0 

MATERIAL 

PLY3 

P100COPR 

0.40 

0.12 

0.01 

00.0 


$ 


Define Design 

Variables 



$ 

$ temperature 

variables 





* 

VARIABLE 

TEMP12 

1 

12 

70.0 

1500.0 

1500 

VARIABLE 

TEMP13 

2 

13 

70.0 

1200.0 

1075 

VARIABLE 

TEMP14 

3 

14 

70.0 

1200.0 

650 

$ link 

temperatures in ply 

1 to 

plies 2 and 3 



DEPEND 

TEMP22 

22 

12 

1.0 



DEPEND 

TEMP23 

23 

13 

1.0 



DEPEND 

TEMP24 

24 

14 

1.0 



DEPEND 

TEMP32 

32 

12 

1.0 



DEPEND 

TEMP33 

33 

13 

1.0 



DEPEND 

TEMP34 

34 

14 * 

1.0 
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$ 


$ ... 
e 

pressure 

variables 





9 

VARIABLE 

PRES1 

4 

111 

-25000. 

0.0 

0. 

VARIABLE 

PRES2 

5 

112 

-25000. 

0.0 

-1000. 

VARIABLE 

PRES3 

6 

113 

-25000. 

0.0 

-1000. 

VARIABLE 

PRES4 

7 

114 

-25000. 

0.0 

-1000. 

§ ensure 

' that N^ 

= Pi * t lf 

throughout tailoring 



DEPEND 

NY1 

201 

111 

0.04 



DEPEND 

NY2 

202 

112 

0.04 



DEPEND 

NY 3 

203 

113 

0.04 



DEPEND 

NY4 

204 

114 

0.04 




$ 

$ interphase properties 


$ DEPEND^ allows all 

plies to 

have equivalent 

interphase characteristics 

9 

VARIABLE modulus 

8 

151 


0.05E08 

3.00E8 0.177E8 

DEPEND 

152 

151 

1.0 



DEPEND 

153 

151 

1.0 



VARIABLE CTE 

9 

161 


0.05E-5 3 

. 00E-5 1.50E-5 

DEPEND 

162 

161 

1.0 



DEPEND 

163 

161 

1.0 



VARIABLE S 

10 

171 


0.05E05 

5.00E5 1.000E5 

DEPEND 

172 

171 

1.0 



DEPEND 

173 

171 

1.0 



VARIABLE td 

11 

141 


0.05 

0.15 0.12 

DEPEND 

142 

141 

1.0 



DEPEND 

143 

141 

1.0 



5“ 1 " - 

S 


Constraints 



5 

$ stress 

S 

$ f iber 

constraints where 

0.85*S C ' < o* < 

0 . 85 * St * 

1.000 

1 

0 

751 

755 

$ longitudinal 
CONSTRNT SfllCl 

&CNST 0.850 

CONSTRNT Sf 11C2 

1.000 

1 

0 

761 

765 

&CNST 0.850 

CONSTRNT Sf 11C3 
&CNST 0.850 

1.000 

1 

0 

771 

775 

$ transverse 
CONSTRNT Sf22Cl 

1.000 

1 

0 

51 

55 

fiCNST 0.850 

CONSTRNT Sf22C2 

1.000 

1 

0 

61 

65 

&CNST 0.850 

CONSTRNT Sf22C3 

1.000 

1 

0 

71 

75 


&CNST 0.850 

S 

$ matrix 

$ longitudinal 

CONSTRNT SmllAl 1.000 1 0 451 455 

&CNST 0.850 

CONSTRNT SmllA2 1.000 1 0 461 465 

fiCNST 0.850 

CONSTRNT SmllA3 1.000 1 0 471 475 

&CNST 0.850 

S 

CONSTRNT SmllBl 1.000 1 0 551 555 

&CNST 0.850 

CONSTRNT SmllB2 1.000 1 0 561 565 

&CNST 0.850 

CONSTRNT SmllB3 1.000 1 0 571 575 

&CNST 0.850 
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$ 

CONSTRNT 

SmllCl 

1.000 

1 

0 

651 

655 

&CNST 

0-850 




661 

665 

CONSTRNT 

SmllC2 

1.000 

1 

0 

&CNST 

0-850 




671 

675 

CONSTRNT 

SmllC3 

1.000 

1 

0 

&CNST 

0.850 






S transverse 




151 

155 

CONSTRNT 

Sm22Al 

1.000 

1 

0 

&CNST 

0.850 




161 

165 

CONSTRNT 

Sm22A2 

1.000 

1 

0 

&CNST 

0.850 




171 

175 

CONSTRNT 

Sm22A3 

1.000 

1 

0 

&CNST 

0.850 






$ 

CONSTRNT 

Sm22B3 

1.000 

1 

0 

251 

255 

&CNST 

0.850 




261 

265 

CONSTRNT 

Sm22B3 

1.000 

1 

0 

&CNST 

0.850 




271 

275 

CONSTRNT 

Sm22B3 

1.000 

1 

0 

&CNST 

0.850 






$ 

CONSTRNT 

Sm22Cl 

1.000 

1 

0 

351 

355 

&CNST 

0.850 




361 

365 

CONSTRNT 

Sm22C2 

1.000 

1 

0 

&CNST 

0.850 




371 

375 

CONSTRNT 

Sm22C3 

1.000 

1 

0 

&CNST 

$ 

$ 

0.850 






interphase 





$ longitudinal 




1051 

1055 

CONSTRNT 

SdllBl 

1.000 

1 

0 

&CNST 

0.850 




1061 

1065 

CONSTRNT 

SdllB2 

1.000 

1 

0 

&CNST 

0.850 




1071 

1075 

CONSTRNT 

SdllB3 

1.000 

1 

0 

&CNST 

0.850 






$ 

CONSTRNT 

SdllCl 

1.000 

1 

0 

1151 

1155 

&CNST 

0.850 




1161 

1165 

CONSTRNT 

SdllC2 

1.000 

1 

0 

&CNST 

0.850 




1171 

1175 

CONSTRNT 

SdllC3 

1.000 

1 

0 

&CNST 

0.850 






S transverse 




851 

855 

CONSTRNT 

Sd22Bl 

1.000 

1 

0 

&CNST 

0.850 




861 

865 

CONSTRNT 

Sd22B2 

1.000 

1 

0 

&CNST 

0.850 




871 

875 

CONSTRNT 

Sd22B3 

1.000 

1 

0 

&CNST 

0.850 






$ 

CONSTRNT 

Sd22Cl 

1.000 

1 

0 

951 

955 

&CNST 

0.850 




961 

965 

CONSTRNT 

Sd22C2 

1.000 

1 

0 

&CNST 

0.850 




971 

975 

CONSTRNT 

Sd22C3 

1.000 

1 

0 

&CNST 

0.850 







$ 

$ lower bounds on 

$ 

CONSTRNT all 

&CNST 0.850 

CONSTRNT a22 


extensional stiffness (85% of initial) 
1.000 2 0 1931 

1.000 2 0 1932 
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&CNST 0.850 

CONSTRNT a33 1.000 2 0 1933 

&CNST 0.850 

$ _ _______ _ 

$ Define optimization analysis control parameters... 


OPTIMIZE 

1 

o 

8 

3440 

1 40 

200 

5 

$ Objective 

Function Data. 





OBJTYPE 

SENSE 

+ 

*END INP 

3 

1.00 1.00 
0.0 

0.00 

0.00 

0.00 0.00 

0.0 0.0+ 
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2.3.2 l±45] s Silicon Carbide/Titanium 


The consolidation temperature and pressure histories are tailored in order to maximize the in- 
service load applied after fabrication for a SiC/Ti 15-3-3-3 composite system. The thermo- 
mechanical (TM) load consists of a linear temperature increase to 600°F followed by the 
application of a biaxial in-plane tensile force (N x = N y = 1500.01b/in) (refer to Figure 2.3-1). 

A 40% FVR is maintained in each ply and a O.Olin ply thickness is assigned to all plies. To 
simulate the fabrication process and subsequent TM loading in the input file (Table 2.3-2), a 
continuation record (+) must be used for the PROCESS, TIME, TEMPRs, PRESR, NX and NY 

records. 


Eight design variables are needed for the simultaneous tailoring of the fabrication process and in- 
service loading: three temperature parameters, four pressure parameters, and the resultant force 
which is maximized. The initial and final temperatures during fabrication are not used along 
with the temperature history used for the in-service loading. The first four control points for the 
pressure are the design variables, the final pressure during fabrication and the in-service pressure 
l oadin g after fabrication are all set to zero. A DEPEND record is used to keep N y equal to e 
design variable N x because only one resultant force is needed as a design variable and the 
DEPEND record ensures equivalency between the loads throughout tailoring. 

Stress constraints are similar to the previous example, except for die elimination of the extra 
microstresses. This is because no interphase is used and the additional stress constraints are 
created by the TM loading after fabrication. If the interphase is non-existent, then only 6 
microstresses exist within the composite model and perfect bonding between the fiber and matrix 
is assumed. Additional constraints are on the final stiffness of the fabricated laminate. 

A different objective function is used for this example, OBJTYPE record identifies the objective 
function to be 4 and the corresponding loads are indicated in the SENSE record. Note that in 
the SENSE record the first two entries correspond to N x and N y , though only N x is a design 
variable, both N x and N y are equalled weighed in the objective function. A -1.0 is assigned m 
the SENSE record indicating maximization of the tensile in-plane force. In contrast, if a 
compressive in-plane force was maximized a positive 1.0 would be used. 
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2.3 Sample Input Data Files 


Table 2.3-2: Input for Case 2 ([±45]s SiC/Til5) 


Material: 

Tailor: 

Objective: 


CASE 2 
MMLT Input Data File 

- [ ±45 ] s Silicon Carbide (SICA) /Titanium (TI15) 

- Fabrication Process 

* Temperature and Pressure Histories 

- Maximize the Tensile Load at an Elevated Temperature 
and Applied After Fabrication 


| Fabrication Process _ 

$~noT”of load history points (NPNT) : no. of subincrements (I(NPNT-l)) 


* 

PROCESS 

8 

10 

20 10 

5 

5 

$ admissible parameters: 

TIME, TEMPR, PRESR, NX, 

NY, NXY , 

$ 

$ time 

(sec) , 

continuation card used 


8000.0 

TIME 


1 0.0 

4000.0 

+ 

13000.0 

14000.0 

18000.0 


a TEMPR 

$ temperature 

history 

(T) of cool-down 

phase , 

TEMPR 



2 1750.0 

1333.0 

866.0 

+ 

300*0 

800.0 

800.0 



TEMPR 


3 1750.0 

1333.0 

866.0 

+ 

300.0 

800.0 

800.0 


866.0 

TEMPR 


4 1750.0 

1333.0 

+ 

300.0 

800.0 

800.0 




10 


10 


400.0 

400.0 

400.0 


70.0+ 

70.0+ 

70.0+ 


0.0 


-40.0 -40.0 


0 . 0 + 

0 . 0 + 

0 . 0 + 


$ consolidation pressure (psi) history 

pj^SH 5 -1000.0 -1000.0 -1000.0 -1000.0 

+ 0.0 0.0 0.0 

$ tensile resultant force (lb/in), in global x-direction 

NX 6 0.0 0.0 0.0 

+ 0.0 0.0 1500.0 

$ resultant force (lb/in), equivalent P ; * t, 

NY 7 -40.0 -40. C 

+ 0.0 0.0 1500.0 __________ 

$ Material and Composite Identification __ 

$ assign ply lay-up 
LAY-UP PLY1 PLY 2 

$ fib/mat 

MATERIAL PLY1 SIC2TI53 

MATERIAL PLY2 SIC2TI53 

MATERIAL PLY 3 SIC2TI53 ================ 

$ Define Design Variables... 

$ 

$ temperature variables 

70.0 1700.0 1333. 

70.0 1700.0 866. 

70.0 1700.0 400. 

1.0 
1.0 


PLY 3 
fvr 

td 

tl 

ang 

0.40 

0.00 

0.01 

45.0 

0.40 

0.00 

0.02 

-45.0 

0.40 

0.00 

0.01 

45.0 


VARIABLE 

TEMP12 

1 

12 

VARIABLE 

TEMP13 

2 

13 

VARIABLE 

TEMP 14 

3 

14 

DEPEND 

TEMP22 

22 

12 

DEPEND 

TEMP23 

23 

13 
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DEPEND 

TEMP24 

24 

14 

1.0 

DEPEND 

TEMP32 

32 

12 

1.0 

DEPEND 

TEMP 3 3 

33 

13 

1.0 

DEPEND 

TEMP 3 4 

34 

14 

1.0 


$ 

$ pressure variables 


5 

VARIABLE 

PRESl 

4 

111 

-50000. 

0.0 

0. 

VARIABLE 

PRES2 

5 

112 

-50000. 

0.0 

-1000. 

VARIABLE 

PRES3 

6 

113 

-50000. 

0.0 

-1000. 

VARIABLE 

PRES4 

7 

114 

-50000. 

0.0 

-1000. 

DEPEND 

NY1 

201 

111 

0.04 



DEPEND 

NY2 

202 

112 

0.04 



DEPEND 

NY 3 

203 

113 

0.04 



DEPEND 

NY4 

204 

114 

0.04 



* 

$ Tensile in- 

i'* 

■plane force 

applied after fabrication at 

an elevated 

VARIABLE 

NX 

8 

198 

500.0 

6500.0 

4000.0 

DEPEND 
$ 

NY 

199 

198 

il 

ii 

II 

II 

II H* 
it * 

II o 

II 

II 

II 

II 

II 

II 

II 

II 

II 

It 

n 

n 

it 

ll 

it 

ll 

n 

n 

ll 

ll 

ll 

n 

ll 

ll 

ll 

ii 

n 

n 

ii 


$ Constraints 

$ stress constraints where 0.85*S c l < a 1 < 0.85*S T l 

$ 

$ fiber 


$ longitudinal 
CONSTRNT SfllCl 

1.000 

1 

0 

751 

758 

&CNST 

0.850 






CONSTRNT 

Sf 11C2 

1.000 

1 

0 

761 

768 

ficCNST 

0.850 






CONSTRNT 

Sf 11C3 

1.000 

1 

0 

771 

778 

&CNST 

0.850 






$ transverse 






CONSTRNT 

Sf 22C1 

1.000 

1 

0 

51 

58 

&CNST 

0.850 






CONSTRNT 

Sf 22C2 

1.000 

1 

0 

61 

68 

&CNST 

0.850 






CONSTRNT 

Sf 22C3 

1.000 

1 

0 

71 

78 

ficCNST 

c 

0.850 






$ 

matrix 






$ longitudinal 
CONSTRNT SmllAl 

1.000 

1 

0 

451 

458 

&CNST 

0.850 






CONSTRNT 

SmllA2 

1.000 

1 

0 

461 

468 

ficCNST 

0.850 






CONSTRNT 

SmllA3 

1.000 

1 

0 

471 

478 

&CNST 

$ 

CONSTRNT 

0.850 






SmllCl 

1.000 

1 

0 

651 

658 

ficCNST 

0.850 






CONSTRNT 

SmllC2 

1.000 

1 

0 

661 

668 

ficCNST 

0.850 






CONSTRNT 

SmllC3 

1.000 

1 

0 

671 

678 

ficCNST 

0.850 






$ transverse 






CONSTRNT 

Sm22Al 

1.000 

1 

0 

151 

158 

ficCNST 

0.850 






CONSTRNT 

Sm22A2 

1.000 

1 

0 

161 

168 

ficCNST 

0.850 






CONSTRNT 

Sm22A3 

1.000 

1 

0 

171 

178 

ficCNST 

0.850 







5 
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2.3 Sample Input Data Files 


CONSTRNT 

&CNST 

CONSTRNT 

&CNST 

CONSTRNT 

&CNST 

Sm22Cl 

0.850 

Sm22C2 

0.850 

Sm22C3 

0.850 

1.000 

1.000 

1.000 

$ 

$ lower 

bounds on 

extensional 

$ 

CONSTRNT 

all 

1.000 

&CNST 

0.850 


CONSTRNT 

a22 

1.000 

&CNST 

0.850 


CONSTRNT 

a33 

1.000 

&CNST 

0.850 


$ ensures proper 

temperature 

CONSTRNT 

DTEMP 

1.000 

&CNST 

2.50 


$ Define optimization analysii 

s, 

OPTIMIZE 

0 

5 

$ Objective Function Data. . . 



OBJTYPE 

4 


SENSE 

-1.00 

-1.00 0. 


o 

0 

0 


351 

358 

361 

368 

371 

378 


0 

0 

0 


1931 

1932 

1933 

:er fc 
1937 


3440 


0*00 


0.00 


1940 


40 


0.00 


200 


0.0 


0.0 


♦END INP 
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PROCEDURE TO RUN MMLT 


3.1 Introduction 

The MMLT computer code has been run in batch mode on the VAX Cluster, Cray X-MP, and 
Convex computer facilities at NASA Lewis Research Center. When using the VAX, the VMS 
operating system is used. In contrast, the Cray X-MP and Convex use UNICOS and UNIX 
operating systems, respectively. The procedure to compile and execute the program on the 
VAX, Cray X-MP, or Convex is given in this section. In each example MMLT. FOR is the 
source code, mmlt.inp is the MMLT input data file, data.dat is the METCAN constituent 
databank, and mmlt or MMLT.EXE is the executable file depending on the system used. Also, 
Table 3.1-1 lists the input/output units needed to read and write to the correct files. The 
compilation and linking of the MMLT code is only required the first time the code is used. 
Afterwards the mmlt or MMLT.EXE can be used every time an execution is desired. 


Table 3.1-1: Read/Write Units for MMLT 


Description 

Input 

Output 

File Name 

MMLT input 

9 


mmlt.inp 

Constituent database 

8 


data.dat 

METCAN output 


16 

metcan.out 

METCAN input 


15 

metcan.inp 

History of ADS Optimization 


7 

h_opt.out 

Postprocessing files (10) 


61-70 

post*. out 

Scratch files (2) 


28-29 

scr*.out 


The * is used as a wildcard symbol. 
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3.2 Run MMLT on the VAX 

The first step is to compile and link the MMLT program to create an executable file named 
MMLT.EXE. The following C_MMLT.COM file can be used to compile the source code: 

! com file to compile and link mmlt .for 
for mmlt 
link mmlt 
exit 

To execute the command file type "C_MMLT" after the VMS prompt ($) on the VAX machine. 
After compiling and linking the code, another COM file has to be created which can be called 
MMLT. COM and is listed below: 


i com 

! 

1 — > input files 
! 

! main input file 
assign mmlt.inp 

! 

1 METCAN database 
assign data.dat 

1 

1 — > output files 
assign metcan.out 
assign metcan.dat 
assign h_opt.out 
J 

1 — > scratch files 
assign scrl.dat 
assign scr2.dat 
! — > post processing 
assign post_l #out 
assign post_2.out 
assign post_3.out 
assign post_4.out 
assign post__5.out 
assign post_6.out 
assign post_7.out 
assign post_8,out 
assign post_9.out 
assign post_10*out 

1 

run NMLT • EXE 


file for MMLT 


for009 


forOOS 


for016 

forOlS 

for007 


for028 

for029 

files from METCAN 
for061 
for062 
for063 
for064 
for065 
for066 
for067 
for068 
for069 
for070 


i 

del/noconf iofl%-dat;* 

del/noconf scr% .dat ; * 

exit 
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3.2 Run MMLT on the VAX 


The following command is issued to start the execution of the code: 

$ SUBMIT MMLT.COM 

where the $ is the VAX prompt. When the execution of the program is completed, the output 
files, as named in the command file, can be reviewed by the user. 
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3.3 Run MMLT on the Cray X-MP 

The job execution file to compile, link, and execute MMLT on the Cray at NASA LeRC can be 
submitted from the VAX or VM front-end machines. Slight differences exist between the two 
files. The procedure used here is very similar to the one used in the previous section (section 
3.1). If using the VAX front-end machine, the file type extension ♦.JOB should be used. The 
listing of MMLT. JOB submitted from the VAX front-end is given as follows: 

# USER=ayuserid PW=»yuserpwd 

# QSUB -r MMLT 

# QSUB -nr 

# QSUB -nc 

# QSUB -q batch 

# QSUB -IT 60 

# QSUB -1M 1.6Mw 

# QSUB -eo 
set -vxk 
cd 

akdir aalt_o 
cd mmlt_o 

j a 

fetch sslt.f -t' melt. for’ 
cft77 -dp -ev melt . f 
segldr -o aalt.exe abit.o 
fetch fort. 09 -t ' mmlt . inp ’ 

fetch fort .08 -t * data . dat 1 

■alt . exe 

dispose fort .16 -t 1 aetcan . out 1 

dispose fort. 15 -t 1 aetcan . inp * 

dispose fort .07 -t * h_opt . out * 

dispose fort. 28 -t' scrl.dat ' 

dispose fort. 29 -fscr2.dat' 

dispose fort. 61 -t *post_l .out • 

dispose fort. 62 — t’post_2.out’ 

dispose fort. 63 -t 1 post_3 . out * 

dispose fort. 64 -t ' post_4 . out 1 

dispose fort. 65 -t ' post_5 . out • 

dispose fort. 66 — t ' post_6 . out ' 

dispose fort. 67 -t ‘post_7 .out 1 

dispose fort. 68 -t 'post_8 .out ’ 

dispose fort. 69 -t 'post_9 .out 1 

dispose fort. 70 -t ’post_10.out • 

ja -s 
exit 
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3.4 Run MMLT on the Convex 


The first step is to compile and link MMLT program to create an executable file named mmlt.o. 
The following mmlt_c file can be used to compile and link the source code: 

echo 

echo "compile? (y/n)" 
read pi 
# 

# compile 

# 

if [ Spl = "y" ] 
then 

echo "COMPILING. . . " 
fc -c -72 -nw -db \ 

mmlt.f mmltin.f metbuf_2.f objf.f cons . f \ 
defn.f auxl.f rept.f mmcv3_l.f u967ads.f 
fi 
# 

# load 

# 

echo "LOADING..." 
fc -db -o mm It \ 

mmlt.o mmltin.o metbuf_2.o mmcv3_l.o objf.o cons.o defn.o auxl.o \ 
rept.o u967ads.o 

# 

echo "END" 

After compiling and linking the code, another executable file should be created which can be 
called mmlt and is listed below: 


# 

# echo 
setenv 
setenv 

# 

# 

# 

setenv 

setenv 

setenv 

# 

# 

# 

setenv 

setenv 

setenv 

setenv 

setenv 

setenv 

setenv 

setenv 

setenv 

setenv 

setenv 

setenv 

# 


input - mmlt 
"SETTING I/O UNITS. . . " 

FOROO 9 mmlt . inp 

FOROO 8 data • dat 

output file(s) 

FORO 16 met can . out 

FORO 15 metcan . inp 

FORO 07 h_opt . out 

for metcan post processing files 

FOR028 scrl.dat 
FOR029 scr2.dat 
FOR061 postal. out 
FOR062 post_2.out 
FOR063 post_3.out 
FOR064 post_4.out 
FOR065 post_5.out 
FOR066 post_6.out 
FORO 6 7 post_7 . out 

FOR068 post_8.out 
FOR069 post_9.out 
FOR070 post_10.out 
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3.4 Run MMLT on the Convex 


§ run executable 

# 

echo "RUNNING, , ," 
wait 
# 

# reset i/o units 

# 

echo "RESETTING I/O UNITS,,," 

unset env FOR005 

unsetenv POR007 

unset env FOR008 

unsetenv FOR019 

unsetenv FOR039 

unsetenv FOR040 

unsetenv FORO 6 * 

unsetenv FOR070 

echo "END" 
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APPENDIX A 


Symbols: 

E 

F(Z) 

G 

I 

i 

M 

N 

P 

Q(Z) 

s 

T 

T m 

T 0 

t 

ta 

t, 

Z 


Nomenclature 


Modulus (psi). 

Objective function. 

Shear modulus (psi). 

Value of load history subincrement, 
i-th value of a constant. 

Moment resultant (lb-in/in). 

Force resultant (lb/in). 

Pressure (psi). 

Inequality constraint. 

Strength (psi). 

Temperature (°F). 

Melting temperature (°F). 

Reference temperature (°F). 

Time (sec). 

Interphase thickness (%). 

Ply thickness (in). 

Design vector. 
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Greek: 

a Coefficient of thermal expansion (in/in/°F). 
6 Strain. 

6 Ply orientation angle (degrees). 

v Poisson’s ratio, 

p Density (lb. in 3 ). 

a Stress (psi). 


Coordinate Systems: 

x,y,z Global (laminate) coordinate system. 
1,2,3 Local (material) coordinate system. 


Subscripts: 

A,B,C 

d 

f 

1 

L 

m 

U 


Constituent microregion identification. 
Interphase. 

Fiber. 

Laminate. 

Lower. 


Matrix 

Upper. 
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Constituent Databank 


MMLT utilizes the same constituent databank used in METCAN. Furthermore, the following 
description in this appendix contains the same information found in Reference 2 but is included here 
for the user’s convenience. The file contains room temperature constituent (fiber, matrix, and 
interphase) properties of various materials. The primary purpose of the databank is to relieve the 
user of the repetitive burden of including this information as part of the primary input file. By 
assigning a unique code name for each constituent in the databank, the user merely has to specify 
the code name in the MATERIAL record of the primary input file to acquire the numerous 
properties for a specific constituent. Also, advantages of using the databank are the flexibility of 
adding other constituents as they become available and updating current constituent materials due 
to current research. The databank is organized into three sections; the first for the fiber properties, 
the second for matrix properties, and the third for the interphase properties. The organization of 
the databank and the properties contained for each constituent is shown in Figure B-l. 

The databank is organized into three sections; each section containing a set of a number of records 
in a predefined order. The first section contains the fiber properties and contains nine records. The 
second section arranges the matrix properties into ten records, while the third section contains the 
interphase properties on nine records. Each record is composed of one or more lines of data whose 
format varies from record to record. Note that the format for the constituent databank records is 
different from the fixed format employed in the primary input file. Each line of data begins with 
a mnemonic. Character, integer, and exponential formats must be right justified, while the floating 
point format can fall anywhere in the appropriate field. The various constituent properties input into 
the databank are room temperature values. The units can be given in English or SI but the user 
must be consistent throughout. 
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Figure B-l: Organization of the Constituent Databank 


Constituent databank 
Contains room temperature material 
properties of the fiber, matrix and interphase 


1 

Fiber material properties 

1 

Matrix material properties 

1 

Interphase material properties 

• Fiber diameter 

• Weight density 

• Weight density 

• Weight density 

• Modulus 

• Interphase thickness 

• Melting temperature 

• Poisson’s ratio 

• Modulus 

• Moduli 

• Coefficient of thermal 

• Poisson’s ratio 

• Poisson’s ratios 

expansion 

• Coefficient of thermal 

• Coefficients of thermal 

• Thermal conductivity 

expansion 

expansion 

• Heat capacity 

• Thermal conductivity 

• Thermal conductivities 

• Strengths 

• Heat capacity 

• Heat capacity 

• Strains 

• Strengths 

• Strengths 

• Void thermal conductivity 

• Strains 

• Reference stresses 

• Melting temperature 

• Void thermal conductivity 

• Reference stress-rates 

• Reference stresses 

• Melting temperature 

• Exponents for multi- 

• Reference stress-rates 

• Reference stresses 

factor interaction 
relationship 

• Exponents for multi- 
factor interaction 
relationship 

• Reference stress-rates 

• Exponents for multi- 
factor interaction 
relationship 
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Appendix B 


B.l Fiber Entries 


B. 1 . 1 Fiber Data Record: CODE 


Description: Defines fiber code name 


Format and Example: 



CODE COMMENT 

+ 

SICA Optional description of fiber 

+ 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

CODE 

A4 

1-4 

Fiber code name 

2 

COMMENT 

A76 

5-80 

Optional description of fiber 


Remarks : 

1. The CODE for each fiber in the databank must be unique. 
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B.1.2 Fiber Data Record: FP 


Description: Defines fiber physical properties 


Format 

and Example: 




FP 

N f 

Df 

Pf 

T M f 


FP 

1 

. 560E-02 

. 110E+00 

•487E+04 


Field 

No. 

Field Name Format 

Columns 

Description 

1 


FP 

A2 

1-2 

Mnemonic name 

2 


N f 

15 

4-9 

Number of fibers per bundle 

3 


D, 

E10.3 

10-19 

Fiber diameter 

4 


Pf 

E10.3 

20-29 

Fiber weight density 

5 


T Mf 

E10.3 

30-39 

Fiber melting temperature 
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B.1.3 Fiber Data Record: FE 


Description: Defines fiber elastic properties 


Format and Example: 


E 

+ — 

E 

+ — 

n 

A 

n 

G 

G 

: . 620E+08 

+- 

. 620E+08 
+ — 

. 300E+00 
+ -• 

. 300E+00 

. 238E+08 

.238E+08 


Field No. 

Field Name 

Format 

Columns; 

Description 

1 

FE 

A2 

1-2 

Mnemonic name 

2 

E 

E10.3 

4-13 

Fiber 

longitudinal 





modulus 

3 

E 

E10.3 

14-23 

Fiber 

transverse modulus 

4 

n 

E10.3 

24-33 

Fiber 

Poisson's ratio 

5 

n 

E10.3 

34-43 

Fiber 

Poisson's ratio 

6 

G 

E10.3 

44-53 

Fiber 

in-plane shear modulus 

7 

G 

E10.3 

54-63 

Fiber 

out-of-plane shear 


modulus 
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B.1.4 Fiber Data Record: FT 


Description: Defines fiber thermal properties 

Format and Example: 

i L — — — — — — — — — — — — — — —— — — 

h — 

FT otm 

— — — — — -r- 

«f22 

Kfll *M22 

<v 


FT . 180E- 

— + — — — — 

■05 . 180E-05 

. 750E+00 . 750E+00 . 

290E+00 

- + 

Field No. 

F +* 

Field Name 

Format 

Columns 

Description 

1 

FT 

A2 

1-2 

Mnemonic name 

2 

a m 

E10.3 

4-13 

Fiber longitudinal 





coefficient of 





thermal expansion 

3 

OCf22 

E10.3 

14-23 

Fiber transverse 





coefficient of 





thermal expansion 

4 

Kfl 1 

E10.3 

24-33 

Fiber longitudinal 





thermal conductivity 

5 

Kf22 

E10. 3 

34-43 

Fiber transverse 





thermal conductivity 

6 

<V 

E10.3 

44-53 

Fiber heat capacity 
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B.1.5 Fiber Data Record: FS 


Description: Defines fiber strength allowable 


Format and Example: 


— F — 


— + 

+ -. 




FS 

s nrr s nic 

S G2T 

Sf22C 

Sf12 

s m 


-+ — 



+ _■ 




FS 

. 500E+06 . 650E+06 

500E+06 

. 650E+06 

. 300E+06 

. 300E+06 



+ — 


+ — 




Field No. 

Field Name 

Format 

Columns 

Description 

1 

FS 

A2 

1-2 

Mnemonic name 

2 

SfllT 

E10.3 

4-13 

Fiber longitudinal 
tensile strength 

3 

Sflic 

E10.3 

14-23 

Fiber longitudinal 
compressive strength 

4 

Sq2T 

E10.3 

24-33 

Fiber transverse 
tensile strength 

5 

Sf22C 

E10.3 

34-43 

Fiber transverse 
compressive strength 

6 


E10.3 

44-53 

Fiber in-plane shear 
strength 

7 

s 

E10.3 

54-63 

Fiber out-of-plane 
shear strength 


1 . m cases where only the experimental tensile strength is available, the 
following approximate values are used for the compressive and shear strengths: 
S c = S T and S s = 0.6 St 
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B.1.6 Fiber Data Record: SIGFO 


Description: Defines reference values of fiber stress 


Format and Example: 

+_ + + + + + 

SIGFO O m Of22 Of\2 a m °tl 3 <*03 

+ + + + + + 

SIGFO 0. 0. 0. 0. 0. 0. 


Field No. 
1 
2 

3 

4 

5 

6 
7 


Field Name 
SIGFO 
<*ni 

C fl2 

C f23 

a fl3 

<*03 


Format 

A5 

E10.3 

E10.3 

E10.3 

E10.3 

E10.3 

E10.3 


Columns 

1-5 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 


Description 
Mnemonic name 
Reference longitudinal 
stress value 
Reference transverse 
stress value 

Reference in-plane shear 
stress value 
Reference out-of-plane 
shear stress value 
Reference out-of-plane 
shear stress value 
Reference through the 
thickness stress value 


Remarks : 

1. The reference values of stress are used in the stress term of the multifactor 
interaction relationship and are typically set to zero for convenience. 
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B.1.7 Fiber Data Record: DOTFO 


Description: Defines reference values of fiber stress rate 


Format and 

Example: 





DOTFO 

a m °m 

°fl2 


a fl3 °D3 


DOTFO 

0 . 0 . 

0 . 

— 

0 . 

0 . 0 . 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

DOTFO 

A5 

1-5 

Mnemonic name 

2 

°ni 

E10.3 

11-20 

Reference longitudinal 
stress rate value 

3 

a f22 

E10.3 

21-30 

Reference transverse 
stress rate value 

4 

a fl2 

E10.3 

31-40 

Reference in-plane shear stress 
rate value 

5 


E10.3 

41-50 

Reference out-of-plane shear 
stress rate value 

6 

a fl3 

E10.3 

51-60 

Reference out-of-plane 
shear stress rate value 

7 

°D3 

E10.3 

61-70 

Reference through the 
thickness stress rate value 

Remarks : 


1. The reference values of stress rates are used in the stress rate term of the 
multifactor interaction relationship and are typically set to zero for 


convenience . 
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B.1.8 Fiber Data Record: EXPONENTS 


Description: Defines values of exponents in the multifactor interaction relationship 


Format and Example: 


+- 

+ - 

+- 

-1 — 

— +- 

+- 

— +- 

+ - 

— +- 

+ - 


EXPONENTS 

FEX1 

FEX2 

FEX3 

FEX4 

FEX5 

FEX6 

FEX7 

FEX8 

FEX9 

FEX10 


+ ~ 


+ - 

+ - 








EXPONENTS 

0* 25 

0.25 

0.25 

0.25 

0.25 

0.00 

0.25 

0.25 

0.25 

0.25 

+ - 




— + . 

+- 

— +- 

— 

— 

— +- 


+- 


+- 

+ . 

+- 

— + . 

— + . 


+- 



EXPONENTS 

FEXll 

FEX12 

FEX13 

FEX14 

FEX15 

FEX16 

FEX17 

FEX18 

FEX19 

FEX20 

+- 

+. 

+ - 

+ - 

+ - 

+ - 



+ - 


+ - 


EXPONENTS 

0.25 

0.25 

0.25 

0.25 

0.00 

0.25 

0.25 

0.00 

0.25 

0.25 

+- 

+ . 

+- 

+ - 

— +- 

+- 

— + . 

+- 



+- 



+. 

+. 

F_ 

— +. 

+- 

— +. 



+- 


EXPONENTS 

FEX21 

FEX22 

FEX23 

FEX24 

FEX25 

FEX26 

FEX27 

FEX28 

FEX29 

FEX30 

+- 

+- 

+ - 

+- 

+. 

+ - 

+- 

+. 


+. 


EXPONENTS 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

+ - 


+. 

+- 

+ - 



+. 

+- 

+. 

+. 


+ - 

.+. 

+- 

+. 

+- 




+- 

+ - 



EXPONENTS 

FEX31 










+- 

+. 

+. 

+. 

+ . 

+. 


+. 


+■ 


EXPONENTS 

0.25 











+. 


+. 

+. 

+. 




+. 
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Field No- 

Field Name 

Format 

Columns 

Line 1 : 




1 

EXPONENTS 

A9 

1-9 

2 

FEX1 

F5.2 

11-15 

3 

FEX2 

F5.2 

17-21 

4 

FEX3 

F5.2 

23-27 

5 

FEX4 

F5.2 

29-33 

6 

FEX5 

F5.2 

35-39 

7 

FEX6 

F5.2 

41-45 

8 

FEX7 

F5.2 

47-51 

9 

FEX8 

F5.2 

53-57 

10 

FEX9 

F5.2 

59-63 

11 

FEX10 

F5.2 

65-69 


Line 2 : 


1 

EXPONENTS 

A9 

1-9 

2 

FEX11 

F5.2 

11-15 

3 

FEX12 

F5.2 

17-21 

4 

FEX13 

F5.2 

23-27 

5 

FEX14 

F5.2 

29-33 

6 

FEX15 

F5.2 

35-39 

7 

FEX16 

F5.2 

41-45 

8 

FEX17 

F5.2 

47-51 

9 

FEX18 

F5.2 

53-57 

10 

FEX19 

F5.2 

59-63 

11 

FEX20 

F5.2 

65-69 


Metal Matrix Laminate Tailoring 
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Description 

Mnemonic name 
Reserved for future use 
Modulus-stress rate exponent 
Modulus-stress exponent 
Modulus-temperature exponent 
Strength-stress rate exponent 
Strength-stress exponent 
Strength-temperature exponent 
Poisson’s ratio-stress rate 
exponent 

Poisson's ratio-stress exponent 
Poisson ' s ratio-temperature 
exponent 

Mnemonic name 

Reserved for future use 

Reserved for future use 

Reserved for future use 

Heat conductivity-stress rate 

exponent 

Heat conductivity-stress exponent 
Heat conductivity-temperature 
exponent 

CTE-stress rate exponent 
CTE-stress rate exponent 
CTE-temperature exponent 
Reserved for future use 
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Line 3: 


1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

FEX21 

F5.2 

11-15 

Reserved 

for 

future 

use 

3 

FEX22 

F5.2 

17-21 

Reserved 

for 

future 

use 

4 

FEX23 

F5.2 

23-27 

Reserved 

for 

future 

use 

5 

FEX24 

F5.2 

29-33 

Reserved 

for 

future 

use 

6 

FEX25 

F5.2 

35-39 

Reserved 

for 

future 

use 

7 

FEX26 

F5.2 

41-45 

Reserved 

for 

future 

use 

8 

FEX27 

F5.2 

47-51 

Reserved 

for 

future 

use 

9 

FEX28 

F5.2 

53-57 

Reserved 

for 

future 

use 

10 

FEX29 

F5.2 

59-63 

Reserved 

for 

future 

use 

11 

FEX30 

F5.2 

65-69 

Reserved 

for 

future 

use 

e 4: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

FEX31 

F5.2 

11-15 

Reserved 

for 

future 

use 
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Remarks : 

1, The various exponents listed are used in the multif actor interaction 

equation (Section 1.6.1). For example FEX4 would be the exponent used to 
determine the temperature effect on the fiber ^ modulus . 


E 


fll 


■"fllo 



FEX4 


2. The values of the exponents listed above represent the default values. 

3. The locations of the currently used exponents are shown in Table B.l-1 for 
conven ient re f er ence . 


Table B.l-1: Location of Exponents in Database File for Fiber 


Exponent Field Name (i.e. FEX ) 

Property 

Stress Rate 

Stress 

Temperature 

Modulus 

2 

3 

4 

Strength 

5 

6 

7 

Poisson* s Ratio 

8 

9 

10 

Heat Conductivity 

14 

15 

16 

Coeff. of Thermal Exp. 

17 

18 

19 


All other locations are reserved for future use. 
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B.1.9 Fiber Data Record: OVER 


Description: Defines the end of the fiber section in the databank 


Format : 


OVER END OF FIBER PROPERTIES 


Field No. Field Name Format Columns Description 

1 OVER... A28 1-28 Mnemonic name 

Remarks : 

1. This record is used to differentiate between the fiber and matrix sections of 
the databank. This record is only inserted after the last fiber entry. 
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B.2 Matrix Entries 


B.2.1 Matrix Data Record: CODE 


Description: Defines matrix code name 


Format and Example: 




CODE COMMENT 



TI15 Optional description of matrix 




Field No. 

Field Name 

Format 

Columns 

Description 

1 

CODE 

A4 

1-4 

Matrix code name 

2 

COMMENT 

A76 

5-80 

Optional description of matrix 


Remarks : 

1. The CODE for each matrix in the databank must be unique. 
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B.2.2 Matrix Data Record: MP 


Description: Defines matrix physical properties 


Format and Example: 


-+ 


+< 


MP p m 


-+ +- 

MP . 172E+00 
-+ +- 


Field No. 
1 
2 


Field Name Format 


MP 

D 


A2 

E10.3 


Columns Description 

1-2 Mnemonic name 

4-13 Matrix weight density 
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B.2.3 Matrix Data Record: ME 


Description: Defines matrix elastic properties 


Format and Example: 


ME E m v m a m 

-+ + + + 

ME . 123E+08 . 320E+00 .450E-05 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

ME 

A2 

M 

1 

to 

Mnemonic name 

2 

E* 

E10.3 

4-13 

Matrix modulus 

3 


E10.3 

14-23 

Matrix Poisson's ratio 

4 

a® 

E10.3 

24-33 

Matrix coefficient of thermal 
expansion 
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B.2.4 Matrix Data Record: MT 


Description: Defines matrix thermal properties 


Format and Example: 

+ +■ 

MT K„ 

-+ + + 

MT . 390E+00 . 120E +00 
-+ + + 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

MT 

A2 

1-2 

Mnemonic name 

2 

K. 

E10.3 

4-13 

Matrix thermal conductivity 

3 

Cpm 

E10.3 

14-23 

Matrix heat capacity 


Metal Matrix Laminate Tailoring 
NASA Lewis Research Center 


Part V - Appendices 

B-18 


Appendix B 


B.2.5 Matrix Data Record: MS 


Description: 

Defines matrix strength properties 

Format and Example: 



MS S mT 

^mC S e mT € znC E mS £ mr 



MS 130E+06 

. 130E+06 . 780E+05 .120E+00 . 120E+00 .120E+00 . 120E+00 


-+ + + + + + + + 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

MS 

A2 

1-2 

Mnemonic name 

2 

s aT 

E10.3 

4-13 

Matrix 

tensile strength 

3 

Smc 

E10.3 

14-23 

Matrix 

compressive strength 

4 

®mS 

E10.3 

24-33 

Matrix 

shear strength 

5 

C mT 

E10.3 

34-43 

Matrix 

tensile strain 

6 

£ mC 

E10.3 

44-53 

Matrix 

compressive strain 

7 

€ m S 

E10.3 

54-63 

Matrix 

shear strain 

8 


E10.3 

64-73 

Matrix 

torsional strain 


Remarks : 

1. In cases where only the experimental tensile strength is available, the 
following approximate values are used for the compressive and shear 
strengths: 

= and S s — 0*6 S-J- 
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B.2.6 Matrix Data Record: MV 


Description: Defines additional matrix thermal properties 


Format and Example: 


-+ — 





MV 

Kv 




-+ 





MV 

.019E+00 

• 180E 

+04 




Field No* 

Field Name 

Format 

Columns 

Description 

1 

MV 

A2 

1-2 

Mnemonic name 

2 


E10.3 

4-13 

Matrix void thermal conductivity 

3 


E10.3 

14-23 

Matrix melting temperature 
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B.2.7 Matrix Data Record: SIGMO 


Description: Defines reference values of matrix stresses 


Format and Example: 



Field No. 
Line 1: 

1 

2 

3 


Field Name 

Format 

Columns 

SIGMO 

A5 

1-5 

a mllA 

E10.3 

11-20 

a m22A 

E10.3 

21-30 


Description 
Mnemonic name 

Reference stress in region mllA 
Reference stress in region m22A 
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4 

^m22B 

E10.3 

31-40 

Reference 

stress 

in 

region 

m 

CN 

CM 

6 

5 

°mZ2C 

E10.3 

41-50 

Reference 

stress 

in 

region 

xn22C 

6 

a mi2A 

E10.3 

51-60 

Reference 

stress 

in 

region 

ml2A 

7 

a ml2B 

E10.3 

61-70 

Reference 

stress 

in 

region 

ral2B 

8 

a ml2C 

E10.3 

71-80 

Reference 

stress 

in 

region 

ml2C 

Line 2: 









1 

SIGMO 

A5 

1-5 

Mnemonic name 




2 

a m23A 

E10.3 

11-20 

Reference 

stress 

in 

region 

m2 3 A 

3 

a m23B 

E10.3 

21-30 

Reference 

stress 

in 

region 

ra23B 

4 

a m23C 

E10.3 

31-40 

Reference 

stress 

in 

region 

m23C 

5 

^mJ3A 

E10.3 

41-50 

Reference 

stress 

in 

region 

ml3A 

6 

a ml3B 

E10.3 

51-60 

Reference 

stress 

in 

region 

m!3B 

7 

a ml3C 

E10.3 

61-70 

Reference 

stress 

in 

region 

ml3C 

8 

°m33A 

E10.3 

71-80 

Reference 

stress 

in 

region 

m33A 

Line 3: 









1 

SIGMO 

A5 

1-5 

Mnemonic name 




2 

a m33B 

E10.3 

11-20 

Reference 

stress 

in 

region 

m33B 

3 

a m33C 

E10.3 

21-30 

Reference 

stress 

in 

region 

m33C 

4 

a milB 

E10.3 

31-40 

Reference 

stress 

in 

region 

mllB 

5 

a mIIC 

E10.3 

41-50 

Reference 

stress 

in 

region 

ml 1C 


Remarks : 

1. The reference values of stress are used in the stress term of the 

multifactor interaction relationship and are typically set to zero for 
convenience . 
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B.2.8 Matrix Data Record: DOTMO 


Description: Defines reference values of matrix stress rates 


Format and Example: 


DOTHO CTmilA °m22A °m22C °ml2A °ml2B °talX 




Field No. 

Field Name 

Format 

Columns 

Description 



Line 1 : 







1 

DOTMO 

A5 

1-5 

Mnemonic name 



2 

a mIlA 

E10.3 

11-20 

Reference stress 

rate 

in region mllA 

3 

°m22A 

E10.3 

21-30 

Reference stress 

rate 

in region m22A 
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4 

c7 m22B 

E10.3 

31-40 

Reference 

m22B 

stress 

rate 

in 

region 

5 

°m l22C 

ElO. 3 

41-50 

Reference 

xn22C 

stress 

rate 

in 

region 

6 

°mJ2A 

E10 • 3 

51-60 

Reference 

m!2A 

stress 

rate 

in 

region 

7 

a ml2B 

ElO. 3 

61-70 

Reference 

ml2B 

stress 

rate 

in 

region 

8 

Q ml2C 

ElO • 3 

71-80 

Reference 

m!2C 

stress 

rate 

in 

region 


Line 2: 


1 

DOTMO 

A5 

1-5 

Mnemonic name 




2 

a m23A 

ElO. 3 

11-20 

Reference 
m2 3 A 

stress 

rate 

in 

region 

3 

a m23B 

ElO. 3 

21-30 

Reference 

m23B 

stress 

rate 

in 

region 

4 

a m23C 

ElO. 3 

31-40 

Reference 

m23C 

stress 

rate 

in 

region 

5 

°mnK 

ElO. 3 

41-50 

Reference 

ml3A 

stress 

rate 

in 

region 

6 

°ml3B 

ElO. 3 

51-60 

Reference 

ml3B 

stress 

rate 

in 

region 

7 

a ml3C 

ElO. 3 

61-70 

Reference 

ml3C 

stress 

rate 

in 

region 

8 

a m33A 

ElO. 3 

71-80 

Reference 
m3 3 A 

stress 

rate 

in 

region 
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Line 3: 


1 

DOTMO 

A5 

1-5 

Mnemonic name 




2 

a m33B 

E10.3 

11-20 

Reference 

m33B 

stress 

rate 

in 

region 

3 

a m33C 

E10.3 

21-30 

Reference 

m33C 

stress 

rate 

in 

region 

4 

a ml]B 

E10.3 

31-40 

Reference 

mllB 

stress 

rate 

in 

region 

5 

a mlJC 

E10.3 

41-50 

Reference 

mile 

stress 

rate 

in 

region 


Remarks : 

1. The reference values of stress rates are used in the stress rate term of 

the multifactor interaction relationship and are typically set to zero for 
convenience. 
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B.2.9 Matrix Data Record: EXPONENTS 

Description: Defines values of exponents for matrix properties in the 

multifactor interaction relationship 

Format and Example: 



+- 

+- 

+- 

+- 

+- 

+- 

+- 

+- 

+- 


EXPONENTS 

ME XI 

MEX2 

MEX3 

MEX4 

MEX5 

MEX6 

MEX7 

MEX8 

MEX9 

MEX10 


+ - 

+- 



+- 

+ . 

— + . 

+ - 

+ - 


EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 

0.50 

0.50 





+ . 

+- 

— + - 

— +- 

+- 

— +- 


+- 

+- 

+- 

— - — +_ 

— +- 

+- 



+- 


+- 


EXPONENTS 

MEXll 

MEX12 

MEX13 

MEX14 

MEX15 

MEX16 

MEX17 

MEX18 

MEX19 

MEX20 

+- 

+- 

+- 

+- 

+- 

+- 

+ - 

+ - 

+ - 

-f- 


EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 

0.00 

0.50 

0.50 

+- 

+- 

+ . 

+- 


+- 

— +- 

— +- 

+- 

+. 


+- 

+- 

+. 

+- 

— +. 

+- 

+- 


+. 

+- 


EXPONENTS 

MEX21 

MEX22 

MEX23 

MEX24 

MEX25 

MEX26 

MEX27 

MEX28 

MEX29 

MEX30 

+- 

+ . 

+ - 

+ . 

+. 

+. 

+- 

+. 

+- 

+. 


EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

+- 



+ - 


+ . 

+- 

+ - 



+ 


— +- 

+ - 

+. 

+ . 


+ . 


+ - 

+. 



EXPONENTS 

MEX31 











+- 


+- 


+ - 

— +• 



+- 



EXPONENTS 

0.50 











+. 

+. 

+. 

+. 


+. 


+. 

+. 
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Field No. 

Field Name 

Format 

Columns 

Line 1: 




1 

EXPONENTS 

A9 

1-9 

2 

MEX1 

F5.2 

11-15 

3 

MEX2 

F5.2 

17-21 

4 

MEX3 

F5.2 

23-27 

5 

MEX4 

F5.2 

29-33 

6 

MEX5 

F5.2 

35-39 

7 

MEX6 

F5.2 

41-45 

8 

MEX7 

F5.2 

47-51 

9 

MEX8 

F5.2 

53-57 

10 

MEX9 

F5.2 

59-63 

11 

MEX10 

F5.2 

65-69 


Line 2: 




1 

EXPONENTS 

A9 

1-9 

2 

MEX11 

F5.2 

11-15 

3 

MEX12 

F5.2 

17-21 

4 

MEX13 

F5.2 

23-27 

5 

HEX 14 

F5.2 

29-33 

6 

MEX15 

F5.2 

35-39 

7 

MEX16 

F5.2 

41-45 

8 

MEX17 

F5.2 

47-51 

9 

MEX18 

F5.2 

53-57 

10 

MEX19 

F5.2 

59-63 

11 

MEX20 

F5.2 

65-69 


Description 

Mnemonic name 
Reserved for future use 
Modulus-stress rate exponent 
Modulus-stress exponent 
Modulus-temperature exponent 
Strength-stress rate exponent 
Strength-stress exponent 
Strength-temperature exponent 
Poisson’s ratio-stress rate 
exponent 

Poisson's ratio-stress exponent 
Poisson ' s ratio-temperature 
exponent 

Mnemonic name 

Reserved for future use 

Reserved for future use 

Reserved for future use 

Heat conductivity-stress rate 

exponent 

Heat conductivity-stress exponent 
Heat conductivity-temperature 
exponent 

CTE-stress rate exponent 
CTE-stress rate exponent 
CTE-temperature exponent 
Reserved for future use 
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Line 3: 


1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

MEX21 

F5.2 

11-15 

Reserved 

for 

future 

use 

3 

MEX22 

F5.2 

17-21 

Reserved 

for 

future 

use 

4 

MEX23 

F5.2 

23-27 

Reserved 

for 

future 

use 

5 

MEX24 

F5.2 

29-33 

Reserved 

for 

future 

use 

6 

MEX25 

F5 .2 

35-39 

Reserved 

for 

future 

use 

7 

MEX26 

F5.2 

41-45 

Reserved 

for 

future 

use 

8 

MEX27 

F5.2 

47-51 

Reserved 

for 

future 

use 

9 

MEX28 

F5.2 

53-57 

Reserved 

for 

future 

use 

10 

MEX29 

F5.2 

59-63 

Reserved 

for 

future 

use 

11 

MEX30 

F5.2 

65-69 

Reserved 

for 

future 

use 

Line 4: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

MEX31 

F5.2 

11-15 

Reserved 

for 

future 

use 
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Remarks : 

1. The various exponents listed are used in the multifactor interaction 

equation* For example MEX3 would be the exponent used to determine the 
stress effect on the matrix modulus* 


"mil 



S m if a mli 


MEX3 


2* The values of the exponents listed above represent the default values* 

3. The locations of the currently used exponents are shown in Table B.2-1 for 
convenient reference. 


Table B.2-1: Location of Exponents in Database File for Matrix 


|| Exponent Field Name (i.e. MEX ) 

Property 

Stress Rate 

Stress 

Temperature 

Modulus 

2 

3 

4 

Strength 

5 

6 

7 

Poisson's Ratio 

8 

9 

10 

Heat Conductivity 

14 

15 

16 

Coeff. of Thermal Exp. 

17 

18 

19 


All other locations are reserved for future use* 
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B.2.10 Matrix Data Record: OVER 


Description: Defines the end of the matrix section in the databank 


Format : 


OVER END OF MATRIX PROPERTIES 


Field No. Field Name Format Columns Description 

1 OVER. . . A28 1-28 Mnemonic name 

Remarks : 

1 .This record is used to differentiate between the matrix and interphase sections 
of the databank. This record is inserted only after the last matrix entry. 
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B.3 Interphase Entries 


B.3.1 Interphase Data Record: (No Mnemonic) 


Description: Defines optional interphase description 


Format and Example : 

+ + + + + + + + + + 

TITLE 

+ + + + + + + + + + 

Description of interphase 

+ + + + + + + + + + 

Field No. Field Name Format Columns Description 

1 TITLE A80 1-80 Interphase description 


Remarks : 

1. Unlike the fiber and matrix, which are distinguished by their code names, 
the interphases are defined by their thickness (see DP record) . 
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B.3.2 Interphase Data Record: DP 


Description: Defines interphase physical properties 


Format and Example: 


-+ + +• 

DP p d t d 

-+ + +■ 

DP . 172E+00 .020E+00 

-+ + +• 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

DP 

A2 

1-2 

Mnemonic name 

2 

Pi 

E10.3 

4-13 

Interphase weight density 

3 

t„ 

E10.3 

14-23 

Interphase thickness 


Remarks : 

1. The input data file uses the interface thickness to select the desired 

interphase. This means that the thickness of the desired interphase must 
be unique. 
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B.3.3 Interphase Data Record: DE 


Description: Defines interphase elastic properties 


Format and Example: 



-+ 

-+ 


■ E d 

v d 

a d 


! * 300E+07 

-+— — — — — - 

• 320E+00 
-+ 

“+ — — 

. 450E-05 
--f 



Field No. 

Field Name 

Format 

Columns 

Description 

1 

DE 

A2 

CM 

i 

rH 

Mnemonic name 

2 

Ed 

E10.3 

4-13 

Interphase modulus 

3 

v d 

E10.3 

14-23 

Interphase poisson's ratio 

4 

«d 

E10.3 

24-33 

Interphase coefficient of thermal 
expansion 
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B.3.4 Interphase Data Record: DT 


Description: Defines interphase thermal properties 


Format and Example: 


-+ 

DT K* 

-+ +- 

Cpd 




_ + 

-+ +- 




DT . 390E+00 

. 120E+00 




_+ 

“ + + - 




Field No- 

Field Name 

Format 

Columns 

Description 

1 

DT 

A2 

1-2 

Mnemonic name 

2 

Kd 

E10-3 

4-13 

Interphase thermal conductivity 

3 

C* 

E10-3 

14-23 

Interphase heat capacity 
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B.3.5 Interphase Data Record: DS 


Description: Defines interphase strength properties 


Format and Example: 


-+ + + + + + + + 


DS S<jy 



S dT 

X X 

€ dC C dS C dr 




T T 


DS . 325E+05 

. 325E+05 

. 195E+05 

. 120E+00 . 

L X 

120E+00 - 120E+00 . 120E+00 




-*r — — — — — — — — 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

DS 

A2 

1-2 

Mnemonic name 

2 

S<n 

E10.3 

4-13 

Interphase tensile strength 

3 

S*c 

E10.3 

14-23 

Interphase compressive strength 

4 

Sds 

E10.3 

24-33 

Interphase shear strength 

5 

€ dT 

E10.3 

34-43 

Interphase tensile strain 

6 

C dC 

E10.3 

44-53 

Interphase compressive strain 

7 

e dS 

E10.3 

54-63 

Interphase shear strain 

8 

6 dr 

E10.3 

64-73 

Interphase torsional strain 


Remarks : 

1. In cases where only the experimental tensile strength is available, the 
following approximate values are used for the compressive and shear 
strengths: 

= Sy and — 0*6 Sy 
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B.3.6 Interphase Data Record: DV 


Description: Defines additional interphase thermal properties 


Format and Example: 


-+ 

DV Kv 

-+ +- 




-+ 

DV . 019E+00 

- + — — + - 

. 180E+04 




“ + 

-H +- 




Field No. 

Field Name 

Format 

Columns 

Description 

1 

DV 

A2 

i-2 

Mnemonic name 

2 

Kv 

E10.3 

4-13 

Interphase void thermal 





conductivity 

3 

^Md 

E10.3 

14-23 

Interphase melting temperature 
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B.3.7 Interphase Data Record: SIGDO 


Description: Defines reference values of interphase stresses 


Format and Example: 


+ 

+- 

+ 

— + 

._+ 



— + 

SIGDO 

°diiB 

Oatjr °d22C 

a dl2B 

a dl2C 

a d23B 

a d23C 


+ 

+- 


— + 

+ 



— + 

SIGDO 

0,0 

0.0 0.0 

0.0 

0.0 


O 

* 

O 

o 

o 



,+ + + + + + + + 


SIGDO 


a dl3B 




a dl3C 

+ 

a d33B a d33C 

a dllC 

+ 

— + 







+ -^ 








SIGDO 


0.0 


0.0 

0.0 0.0 

0.0 
















Field 

No. 

Field 

Name 

Format 

Columns 

Description 




Line ! 

I: 











1 

SIGDO 


A5 

1-5 

Mnemonic name 





2 

a dllB 


E10.3 

11-20 

Reference 

stress 

in 

region 

dllB 


3 

a d22B 


E10.3 

21-30 

Reference 

stress 

in 

region 

d22B 


4 

°d22C 


E10.3 

31-40 

Reference 

stress 

in 

region 

d22C 


5 

a dl2B 


E10.3 

41-50 

Reference 

stress 

in 

region 

dl2B 


6 

a dI2C 


E10.3 

51-60 

Reference 

stress 

in 

region 

dl2C 


7 

a d23B 


E10.3 

61-70 

Reference 

stress 

in 

region 

d23B 


8 

°d23C 


E10.3 

71-80 

Reference 

stress 

in 

region 

d23C 
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Line 2 : 


1 

SIGDO 

A5 

1-5 

Mnemonic name 




2 

a dl3B 

E10.3 

11-20 

Reference 

stress 

in 

region 

dl3B 

3 

a dl3C 

E10.3 

21-30 

Reference 

stress 

in 

region 

dl3C 

4 

a d33B 

E10.3 

31-40 

Reference 

stress 

in 

region 

d33B 

5 

a d33C 

E10.3 

41-50 

Reference 

stress 

in 

region 

d33C 

6 

a diic 

E10.3 

51-60 

Reference 

stress 

in 

region 

dllC 


Remarks : 

1. The reference values of the stresses are used in the stress term of the 
multif actor interaction relationship and are typically set to zero for 
convenience. 
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B.3.8 Interphase Data Record: DOTDO 


Description: Defines reference values of interphase stress rates 


Format and Example: 



Field No. Field Name Format Columns Description 

Line 1: 

1 DOTDO A5 1-5 Mnemonic name 

2 a d11B E10.3 11-20 Reference stress rate in region 

dllB 

3 Qxry* E10.3 21-30 Reference stress rate in region 

d22B 

4 °& 2 c E10.3 31-40 Reference stress rate in region 

d22C 

5 o dl2B E10.3 41-50 Reference stress rate in region 

dl2B 
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6 

a dl2C 

E10.3 

51-60 

Reference 

stress 

rate 

in 

region 





dl2C 





7 

Q<\r\\K 

E10.3 

61-70 

Reference 

stress 

rate 

in 

region 





d23B 





8 

a d23C 

E10.3 

71-80 

Reference 

stress 

rate 

in 

region 





d23C 





Line 2 * 
1 
2 

DOTDO 

a dl3B 

A5 

E10.3 

I- 5 

II- 20 

Mnemonic name 
Reference stress 

rate 

in 

region 





dl3B 





3 

a dl3C 

E10.3 

21-30 

Reference 

stress 

rate 

in 

region 





d!3C 





4 

a d33B 

E10.3 

31-40 

Reference 

d33B 

stress 

rate 

in 

region 

5 

a d33C 

E10.3 

41-50 

Reference 

stress 

rate 

in 

region 





d33C 





6 

a d11C 

E10.3 

51-60 

Reference 

dllC 

stress 

rate 

in 

region 


1 . The reference values of the stress rates are used in the stress rate term 
of the multifactor interaction relationship and are typically set to zero 
for convenience. 
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B.3.9 Interphase Data Record: EXPONENTS 

Description: Defines values of exponents for interphase properties in 

multifactor interaction relationship 


+. 

+- 

+. 

+- 

+- 

+- 

+. 

+ - 

+- 

+. 


+ 

EXPONENTS 

DEX1 

DEX2 

DEX3 

DEX4 

DEX5 

DEX6 

DEX7 

DEX8 

DEX9 

DEX10 


+. 

+. 

+ - 

+ - 

+- 

+- 

+. 

+- 

+ - 

+ . 

+ 

+ 

EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 

0.50 

0.50 


+. 

+- 

+- 

+- 

+- 

+- 

+- 


+- 

— + . 

+ 

+ 

+- 

+ . 

+- 

+- 

+- 

+- 

+- 

+- 

+- 

+- 

+ 

+ 

EXPONENTS 

DEXll 

DEX12 

DEX13 

DEX14 

DEX15 

DEX16 

DEX17 

DEX18 

DEX19 

DEX20 


+- 

+ - 

+ - 

+ - 

+- 

+- 

+ - 

+- 

+ - 

+ - 

+ 

+ 

EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 

0.00 

0.50 

0.50 


+- 

+- 

+- 

— + . 

+. 

+- 



+- 

+- 

+- 

+ 

+ 

+- 


— +. 

+- 

+. 

+. 

+- 

— +. 

+. 


+ 

+ 

EXPONENTS 

DEX21 

DEX22 

DEX23 

DEX24 

DEX25 

DEX26 

DEX27 

DEX28 

DEX29 

DEX30 


+- 

+_ 



+. 

+. 

+ . 

+ - 

+ - 

+ - 

+. 

+ 

+ 

EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 


+- 

+ - 

+ - 


+. 


+- 


+ - 

+ - 

+ 

+ 

. +- 

+ - 

+ - 

+- 

+ - 

+ - 

+ - 



+. 

+ - 

+ 

+ 

EXPONENTS 

DEX31 











+- 

+- 


— + . 

+- 



+ . 

+- 

+ - 


+ 

EXPONENTS 

0.50 











+- 

+- 

+. 

+- 

— +- 

+- 

+ - 

+- 

+- 

+ . 
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Field No* 

Field Name 

Format 

Columns 

Line 1: 




1 

EXPONENTS 

A9 

1-9 

2 

DEX1 

F5.2 

11-15 

3 

DEX2 

F5.2 

17-21 

4 

DEX3 

F5.2 

23-27 

5 

DEX4 

F5.2 

29-33 

6 

DEX5 

F5.2 

35-39 

7 

DEX6 

F5.2 

41-45 

8 

DEX7 

F5.2 

47-51 

9 

DEX8 

F5.2 

53-57 

10 

DEX9 

F5.2 

59-63 

11 

DEX10 

F5-2 

65-69 


Line 2: 




1 

EXPONENTS 

A9 

1-9 

2 

DEX11 

F5.2 

11-15 

3 

DEX12 

F5.2 

17-21 

4 

DEX13 

F5.2 

23-27 

5 

DEX14 

F5.2 

29-33 

6 

DEX15 

F5.2 

35-39 

7 

DEX16 

F5.2 

41-45 

8 

DEX17 

F5.2 

47-51 

9 

DEX18 

F5.2 

53-57 

10 

DEX19 

F5.2 

59-63 

11 

DEX20 

F5.2 

65-69 


Description 

Mnemonic name 
Reserved for future use 
Modulus-stress rate exponent 
Modulus-stress exponent 
Modulus-temperature exponent 
Strength-stress rate exponent 
Strength-stress exponent 
Strength-temperature exponent 
Poisson's ratio-stress rate 
exponent 

Poisson's ratio-stress exponent 
Poisson ’ s ratio-temperature 
exponent 

Mnemonic name 

Reserved for future use 

Reserved for future use 

Reserved for future use 

Heat conductivity-stress rate 

exponent 

Heat conductivity-stress exponent 
Heat conduct ivity-temperature 
exponent 

CTE-stress rate exponent 
CTE-stress rate exponent 
CTE-temperature exponent 
Reserved for future use 
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Line 3: 


1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

DEX21 

F5.2 

11-15 

Reserved 

for 

future 

use 

3 

DEX22 

F5.2 

17-21 

Reserved 

for 

future 

use 

4 

DEX23 

F5.2 

23-27 

Reserved 

for 

future 

use 

5 

DEX24 

F5.2 

29-33 

Reserved 

for 

future 

use 

6 

DEX25 

F5 . 2 

35-39 

Reserved 

for 

future 

use 

7 

DEX26 

F5 . 2 

41-45 

Reserved 

for 

future 

use 

8 

DEX27 

F5.2 

47-51 

Reserved 

for 

future 

use 

9 

DEX28 

F5 . 2 

53-57 

Reserved 

for 

future 

use 

10 

DEX29 

F5.2 

59-63 

Reserved 

for 

future 

use 

11 

DEX30 

F5.2 

65-69 

Reserved 

for 

future 

use 

Line 4: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

DEX31 

F5.2 

11-15 

Reserved 

for 

future 

use 
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Remarks : 

1 . The various exponents listed are used in the multifactor interaction 

equation. For example DEX4 would be the exponent used to determine the 
stress effect on the matrix modulus. 


1 

co* 

1 

DEX3 

1 J*dlf a dll. 



2. The values of the exponents listed above represent the default values. 

3. The locations of the currently used exponents are shown in Table B.3-1 for 
convenient reference. 


Table B.3-1: Location of Exponents in Database File for Interphase 


| Exponent Field Name 

(i.e. DEX ) 


Property 

Stress Rate 

Stress 

Temperature 

I Modulus 

2 

3 

4 

I Strength 

5 

6 

7 

Poisson’s Ratio 

8 

9 

10 

Heat Conductivity 

14 

15 

16 

1 Coeff. of Thermal Exp. 

17 


18 

19 


All other locations are reserved for future use 
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